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ABSTRACT 

In  direct-sequence  spread-spectrum  systems,  successful  communications  require  phase 
synchronization  of  the  incoming  pseudonoise  (PN)  coded  waveform  with  a  locally  generated 
replica  at  the  receiver.  It  has  been  previously  shown  that  sequential  PN  code  acquisition  schemes 
have  the  potential  to  achieve  the  best  performance,  but  they  are  tl^  least  analyzed  because  of 


the  analytical  difficulties. 

The  acquisition  time  for  a  PN  code  acquisition  scheme  is  an  important  parameter  for  system 
design  puiposes.  This  thesis  investigates  the  peiformance  of  two  acquisition  schemes  in  terms 
of  the  acquisition  time.  A  fixed  sample  size  (FSS)  test  and  a  truncated  sequential  probability 
ratio  test  (TSPRT)  are  studied  with  noncoherent  demodulation  in  a  classical  additive  white 
Gaussian  noise  (AWGN)  channel  and  in  presence  of  fading.  Optimal  selection  of  desired 
detection  and  false  alarm  probabilities,  the  effects  of  penalty  time,  majority  logic  verification 
schemes  and  channel  signal  to  noise  ratio  (SNR)  mismatch  problems  are  thoroughly  studied  usii^ 


the  flow  graph  technique.  Our  results  show  that  the  TSPRT  is  efficient,  robust  (against  fading), 

fast  and  suitable  for  real  time  low  cost  implementations,  f^cesion  Fy _ .  -- 
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I.  INTRODUCTION 


A.  GENERAL 

The  field  of  spread-spectrum  communications  has  been  around  for  several  decades. 
Although  almost  all  the  applications  prior  to  the  80’ s  fell  in  the  military  domain,  recent 
band  allocation  policies  [Ref.  1]  and  commercially  available  IC  components  have  made 
the  area  appealing  to  commercial  2^)plications  as  well. 

The  basic  signal  characteristics  of  modem  spread-spectrum  systems  can  be  defined 
as  follows  [Ref.  2]. 

1.  The  carrier  is  an  unpredictable,  or  pseudorandom,  wideband  signal. 

2.  The  bandwidth  of  the  carrier  is  much  wider  than  the  bandwidth  of  the  data 

modulation. 

3.  Reception  is  accomplished  by  cross  correlation  of  the  received  wideband  signal 

with  a  synchronously  generated  replica  of  the  wideband  carrier. 

Spread-spectrum  signals  provide  the  following  performance  attributes  [Ref.  3]. 

1.  Low  density  power- spectra  for  signal  hiding. 

2.  Interference  rejection  and  anti  jam  properties. 

3.  Selective  addressing  and  code  division  multiple  access  capability. 

4.  Message  privacy. 

5.  High  resolution  and  ranging. 
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Spread-spectnim  systems  can  be  classified  into  [Ref.  4]  "direct  sequrace", 
"frequency  hopping",  "time  hopping",  "chiip"  and  "hybrid  systems"  with  respect  to  the 
modulation  techniques  used  in  the  system. 

In  direct  sequence  spread-spectrum  systems  (DS-SS),  spectrum  spreading  is 
accomplished  by  modulating  a  data  modulated  signal  a  second  time  using  a  very 
wideband  spreading  signal.  This  second  modulating  wideband  signal  is  chosen  to  have 
properties  which  facilitate  demodulation  of  the  transmitted  signal  by  the  intended  receiver 
and  make  demodulation  by  an  unintended  receiver  as  difficult  as  possible.  This  is  best 
accomplished  by  using  a  signal  that  appears  random  to  the  unintended  receivers  and  that 
can  be  reproduced  by  deterministic  means  in  the  intended  receiver.  Therefore,  this 
waveform  is  usually  referred  as  pseudorandom  noise  (PN)  spreading  signal. 

In  a  spread-spectnim  system  using  frequency  hopping  (FH),  the  carrier  frequency 
is  varied  pseudorandomly  with  time,  whereas  a  time  hopping  (TH)  system  uses 
pseudorandom  time  slots  to  transmit  the  signal.  Hybrid  systems  consist  of  a  combination 
of  two  or  all  three  of  the  DS,  FH,  TH  systems.  The  chirp  spread  spectrum  method  can 
be  considered  as  an  analog  direct  spreading  method  that  does  not  use  a  PN  code. 

B.  PN  CODE  SYNCHRONIZATION 

Code  synchronization  is  vital  for  most  spread-spectrum  systems,  since  the  PN  code 
is  used  as  a  key  for  spreading  and  despreading  the  desired  information.  Throughout  the 
development  of  spread-spectrum  systems,  it  is  true  that  [Ref.  3;  pp.  2141; 
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More  time,  effort  and  money  have  been  spent  developing  and  improving 
synchronizing  techniques  than  in  any  other  area  of  ^read-^)ectrum  systems.  There 
is  no  reason  to  suspect  that  this  will  continue  to  be  true  in  the  future. 

A  complete  coverage  of  the  synchronization  of  spread-spectrum  systems  can  be 
found  in  [Ref.  5].  In  this  thesis,  we  restrict  ourselves  to  direct  sequence  ^read-spectrum 
systems. 

The  data  is  recovered  from  a  DS-SS  communication  system  by  neutralizing  the 
effects  of  the  PN  sequence  that  is  used  for  spreading  the  signal  at  the  transmitter’s  site. 
This  process  is  termed  as  despreading.  Removing  the  PN  sequence  from  the  incoming 
signal  is  primarily  accomplished  by  multiplying  the  incoming  signal  by  a  locally 
generated  and  phase  synchronized  replica  of  the  incoming  PN  sequence.  Therefore,  a 
primary  function  of  the  receiver  is  to  align  the  phase  of  the  local  replica  with  the 
incoming  signal’s  PN  sequence.  The  phase  is  determined  by  correlating  the  two  PN 
sequences.  The  alignment  process  is  usually  accomplished  in  two  stages.  The  first  stage, 
called  acquisition,  brings  the  PN  sequences  into  coarse  alignment,  v’hile  the  second 
stage,  called  tracking,  brings  the  PN  sequences  into  precise  alignment. 

We  study  noncoherent  acquisition  of  the  direct  sequence  spread-spectrum  signals. 
In  this  context  noncoherent  refers  to  the  modulating  carrier,  not  the  PN  sequence. 
Although  coherent  systems  are  simpler  and  more  widely  studied,  they  are  not  practical 
because  code  sequence  acquisition  is  usually  performed  before  recovering  the  modulating 
carrier. 

Since  the  acquisition  process  involves  searching  through  the  uncertainty  phases  of 
the  PN  sequence,  acquisition  schemes  can  be  classified  into  parallel,  serial,  and  hybrid 
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schemes  with  respect  to  search  strategies.  The  parallel  scheme  is  the  fastest  search 
strategy.  A  parallel  scheme  inspects  all  the  uncertainty  phases  simultaneously  and  decides 
which  is  the  most  likely  one.  Each  uncertainty  phase  is  investigated  by  a  path  consisting 
of  a  correlator  and  a  matched  filter.  If  the  period  of  the  PN  sequence  is  large,  the 
hardware  requirement  for  such  a  system  is  excessive  and  hence  is  impractical.  On  the 
other  hand,  a  serial  search  scheme  inspects  one  uncertainty  phase  at  a  time  and 
determines  whether  the  PN  sequences  are  in  alignment.  Hardware  requirements  of  such 
schemes  are  small,  but  the  complexity  of  the  control  process  is  high  and  the  decision 
process  requires  a  longer  time.  Various  combinations  of  the  serial  and  parallel  search 
strategies  are  also  possible. 

A  serial  search  acquisition  detection  scheme  that  uses  a  single  detector  to  examine 
each  of  the  possible  waveform  alignments  for  a  fixed  period  of  time  until  the  correct  one 
has  been  located,  is  termed  as  a  single  dwell  scheme.  If  the  scheme  uses  cascaded 
multiple  detectors  each  having  a  longer  examination  period  than  its  predecessor  to  verify 
the  synchronization  condition,  it  is  called  as  a  multiple  dwell  scheme.  The  test  terminates 
if  any  of  the  detectors  decides  a  non-synchronization  condition.  The  synchronization 
condition  is  accepted  if  all  the  detectors  decide  favorably.  In  this  case  the  examination 
time  is  increased  with  discrete  steps,  and  since  most  of  the  positions  will  be  due  to  a  non¬ 
synchronization  condition,  it  is  hoped  that  drey  wiU  be  detected  and  rejected  in  earlier 
stages  using  a  shorter  time.  These  schemes  both  faU  into  the  category  of  fixed  dwell  time 
or  fixed  sample  size  (FSS)  schemes.  Sequential  schemes  use  a  single  detector  with 
mult^le  thresholds,  therefore  yielding  a  varudrle  examination  time.  It  has  bera  shown  [Ref.  6] 
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that  for  a  two-hypothesis  test  of  a  sampled  random  variable,  the  sequential  probability 
ratio  test  (SPRT)  is  optimum  in  the  sense  that  for  a  given  false  alarm  (Pj^ )  and  d^ection 
(Pj)  probability,  it  requires  minimum  average  number  of  samples  to  produce  a  decision 
if  the  samples  of  the  random  variable  are  independent  and  identically  distributed  (i.i.d.)- 
Therefore,  substantial  savings  can  be  achieved  on  the  overall  acquisition  time  using 
comparable  hardware  to  realize  the  system  where  low-cost,  fast  and  reliable  acquisition 
schemes  are  a  necessity,  especially  in  mobile  af^lications  [Ref.  1], 

This  thesis  is  based  on  the  theory  developed  in  [Ref.  7]  and  [Ref.  8],  two  papers 
which  studied  the  noncoherent  sequential  acquisition  of  PN  sequences  in  terms  of  single 
search  state  basis,  and  the  mean  acquisition  time  respectively.  These  pa^rs  served  as  the 
original  motivation  and  starting  point  for  this  research.  In  this  thesis,  we  also  investigate 
the  variance  of  the  acquisition  time,  and  perform  a  more  detailed  acquisition  time 
analysis  on  the  system  parameters. 

The  rest  of  the  thesis  is  organized  as  follows.  Chapter  n  presents  the  PN  code 
acquisition  system  model  and  derives  the  equations  to  calculate  the  mean  and  the  variance 
of  the  acquisition  time  using  the  flow  graph  technique.  Chapter  m  is  devoted  to  the 
acquisition  time  analysis  on  a  single  dwell  FSS  scheme  and  a  truncated  sequential 
probability  ratio  test  (TSPRT)  scheme  over  some  important  design  considerations. 
Chapter  IV  examines  the  performance  of  both  FSS  and  TSPRT  schemes  under  slowly 
varying  Ricean  fading  conditions.  Finally,  the  conclusions  are  given  in  Chapter  V. 
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n.  SYSTEM  MODEL  AND  FLOW  GRATOS 


In  this  chapter,  we  present  the  acquisition  receiver  system  model*  and  derive  the 
equations  to  determine  the  statistics  of  the  acquisition  time  using  the  flow-graph  technique 
developed  in  [Ref.  9]. 

A.  SYSTEM  MODEL 

Figure  2. 1  depicts  a  block  diagram  of  the  receiver’s  acquisition  system.  We  assume 
that  there  is  no  data  modulation  during  the  acquisition  process.  The  chaimel  is  also 
assumed  to  have  additive  white  Gaussian  noise  with  two-sided  power  spectral  density 
N(/2.  The  input  signal  at  the  receiver  is 

tif)  =  +  6)  +  nit)  (2.1) 

where  y4o  is  the  signal  amplitude,  aft)  is  the  m-sequencv  signal  waveform  with  phase  iAT^ 
(/  is  taken  to  be  an  integer  without  loss  of  generality),  is  the  chip  duration,  A  is  the 
value  determining  how  much  the  timing  of  the  local  PN  generator  is  updated  during  the 
acquisition  process,  and  6  are  the  frequency  and  the  phase  of  the  carrier  and  n(t)  is 
the  additive  noise.  The  local  replica  of  the  m-sequence  has  the  form  aft  +  (/  +  y)ATc), 
where  j  is  an  integer  and  |  7  |  ^  O.S.  The  received  signal  rft)  is  first  despread  by 


'  This  chqiter,  which  presents  the  model  and  framewoA  to  be  used  in  the  subsequoit  clu4>ters,  is  based  on  the 
theoiy  developed  in  [Ref.  7]  and  [Ref.  8J.  The  system  model,  equations  and  design  approximations,  the  flow  graph 
of  the  system  model  and  some  of  the  relevant  equations  derived  from  the  gnqjhs,  and  the  early  termination  feature 
of  the  coincidence  detector  were  originally  propounded  in  these  references.  This  material  is  adapted  and  used  with 
the  permission  of  the  authors. 
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Figure  2.1  Block  diagram  of  noncoherent  serial  acquisition  scheme. 

multiplying  with  the  local  PN  waveform  and  is  then  noncoherently  demodulated.  The  test 
statistics  Y„  is  used  by  the  decision  processor  to  test  if  the  local  and  the  incoming  PN 
sequences  are  aligned  to  within  one  chip  duration.  If  not,  the  local  PN  sequence  phase 
is  updated  by  seconds  and  the  process  repeats.  If  coarse  alignment  is  achieved  then 
i  must  be  equal  to  j  and  tracking  circuitry  takes  over  to  reduce  y  to  zero.  Using  the 
notation  in  Figure  2.1,  after  integrating  over  n  chips  and  ignoring  the  double  frequency 
terms,  the  in  phase  and  quadrature  components  are  given  by 
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(2.2) 


|r(/)a(/+(/n)Ar^cos(u)or)dr  =  ^r^„cosd  +  N,^, 

A 

|r(/)fl(/+(/n)Ar^sin(o,^)dr  =  -^TJS„smd  + 


re^)ectively,  where 


4 

"  |/i(/)a(/+(/+7)A7>»s(Wor)<*, 

"n 

=  |n(/)a(r+(/-Hy)Ar^)sin(wof)dr 


are  independent  Gaussian  random  variables  with  zero  mean  and  variance  =  nT^^A. 


The  expression  for  is 


»r. 

s,  -  lf<i(/.(Arx/H/n)Ar,)  <*. 


For  a  fixed  n,  and  are  indq)endent  Gaussian  random  variables  with 
variance  a„*  and  means  {A(/2)T^„gos,B  and  {A,/2)T^„%vaB  respectively.  Note,  however, 
that  X,  „  and  Xf^  are  not  indq)endent  for  n  ^  m,  and  similarly  for  X,  „  and  X^_„.  The  test 
statistics  for  deciding  alignment  or  non  alignment  is 


The  random  variable  Y„  has  a  non-central  Chi-squared  probability  density  function 
(pdf)  with  two  degrees  of  freedom,  given  by 
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AW  =  y.20 

2tr:  al 


(2.6) 


where  X„  =  [(Ao/2)TJS„cosff]^  +  and  VJ  is  the  modified 

Bessel  function  of  order  zero. 

Based  on  the  test  statistics  Y^,  the  decision  processor  must  decide  wh^her  the  two 
phases  are  aligned  to  within  ^TJ2  of  each  (^r,i.e.,  y  =  /,  or  they  differ  by  at  least  one 
chip.i.e.,  |(/+7)  -  /|  AT^  ^  T^.  For  simplicity,  we  let  i  =  0.  Equivalently,  the  decision 
processor  has  a  task  of  testing  the  following  hypotheses 


{non-alignment)  :  |y+7|  ^  and  jy]  ^1, 
//,  (alignment)  :  y=i=0  and 


(2.7) 


Since  |  7  |  1/2,  under  the  hypothesis  Hq  we  must  have  |  j  \  S:(l/A  -  1/2). 

Note  that  it  is  possible  to  have  (1/2)  <  |  j+7  |  <  (1/A),  which  corresponds  to  neither 
Ho  nor  Hj.  It  falls  between  Hg  and  Hj. 

The  parameter  \  takes  on  different  values  under  Hq  and  Hj.  We  designate 
X„ ,1  as  the  worst  case  values  of  X„  under  Hg  and  H,  respectively.  Using  these  worst-case 
parameter  values,  the  likelihood  ratio  for  (2.7)  can  be  written  as 


A.W 


AC.Iw.) 

frp.W 


o.» 


We  will  now  derive  the  nominal  worst  case  values  for  Xg,  assuming  an  m-sequence 
is  used  as  the  PN  code  sequence  of  communication  system  in  consideration.  Denote  the 


chips  of  the  PN  sequence  by  c*,  where  k  =  ...,-1,0, 1,2,...,  and  c*  =  ±1.  Lrt  iV=2“- 
1  be  the  period  of  the  m-sequence.  The  PN  waveform  is 


a(t)  =  5^  cj)^{t-kT^  (2.9) 

where  PtcU)  is  a  rectangular  pulse  of  amplitude  1  and  width  T^.  Assuming  i  =  0,  we 
obtain  an  expression  for  S„  from  (2.4)  and  (2.9)  as 


Ec*[(l -  Iy  I A)c*  +  (Iy  I A)c*,^  )] 
*-o 

k-O 

/i(1-|y|A)  +  (|y|A)*:Cc^*,^y) 

d-l  »-l 

*•0  k-C 


under 


under 


(2.10) 


where  sgn(x)  is  1  for Jt  2:  0  and  -1  for jc  <  0,  /=  L0+y)^J  6=(j+7)A-  L(3+y)^J 

Note  that  I  or  -1  and  0^  5  <1. 

Deflning  the  per-chip  signal  to  noise  ratio  as 


SNR  = 


2A/o’ 


(2.11) 


we  have 
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(2.12) 


Usrnsli  under  H. 
n 

Usmslo  under  H. 
n 

Since  the  exact  values  of  5„  ^  and  5,  /  are  unknown  to  the  receiver  in  advance,  some 
nominal  worst  case  values  must  be  used  in  designing  the  system.  Results  from 

simulations  for  m-sequences  suggests  that  modeling  ••  0  under  Hi,  and 


modeling  and  V<^*<^**i*;  ^  under //q  yield  the  desired  results.  Using  these 

lk*0 

approximations  we  have 


2<r 


»  n(SiVR)(l-|7|A)^ 
i  SNR 


under  H. 

lo\ 

under 

2a\ 


(2.13) 


These  nominal  worst  case  values  will  be  used  in  the  design  methods  of  various 
acquisition  schemes.  The  design  methods  for  fixed  dwell  and  sequential  schemes  will  be 
discussed  in  later  chapters. 

Regardless  of  the  acquisition  scheme,  some  type  of  verification  logic  is  usually 
incorporated  in  these  systems  to  counter  balance  the  effects  of  costly  false  alarms.  Here, 
we  employ  a  majority  logic  verification  scheme  such  that  in  case  of  a  hit,  (deciding 
is  true)  A  or  less  additional  tests  are  performed  before  committing  a  final  decision.  The 


initial  decision  will  hold  if  B  i^A)  additional  tests  are  favorable;  otherwise,  it  is 
overturned  and  the  search  continues  for  the  next  searching  state.  We  also  assume  that  a 
counter  is  available  and  the  coincidence  detector  (CD)  employs  an  early  termination 
mechanism  in  order  to  save  time.  Particularly,  the  test  will  be  terminated  once  B 
favorable  tests  are  accumulated  or  when  there  is  no  chance  that  B  favorable  tests  can  be 
obtained.  Figure  2.2  describes  the  coincidence  detection  process  with  early  termination 
in  detail. 

B.  FLOW  GRAPH  OF  THE  ACQUISITION  SCHEME 

In  [Ref.  9],  a  flow-graph  technique  was  proposed  for  determining  the  statistics  of 
the  acquisition  time  for  serial  search  acquisition  schemes.  A  circular  flow-graph  was 
developed  using  the  mailcovian  nature  of  the  underlying  serial  search  acquisition  process. 
We  will  briefly  discuss  the  technique  in  the  following. 

The  basic  idea  behind  this  technique  comes  from  deHning  each  possible  relative 
position  of  the  phases  of  the  two  PN  sequences  (incoming  and  local  rqilica)  as  a  state  of 
a  discrete  maricov  process.  There  will  be  a  finite  number  of  states  dqiending  the  value 
of  the  A  (the  parameter  which  determines  how  much  the  phase  of  the  local  PN  sequence 
is  advanced  after  each  Hq  decision). 

Let  Py(n)  indicate  the  probability  that  the  Mailcov  process  will  move  from  state  i 
to  state j  inn  steps  and  let  z denote  the  unit  delay  operator.  The  state  transition  diagram 
can  be  mapped  into  its  equivalent  flow  graph  if  each  transition  branch  from  i  to  j  in  the 
Markovian  diagram  is  assigned  a  gain  equal  to  where  p,  A  p^(l)  is  the  one-step 
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INITIAL  DECISION 


Figure  2.2  Coincidence  detection  with  early  termination. 

transition  probability  and  z  represents  the  unit  delay  associated  with  that  transition.  The 
generating  function  for  this  transformation  is 


P^iz)  =^p^in)z’'.  <2.14) 

Note  that  if  the  time  delay  associated  with  the  transitions  is  specifically  t  seconds 
z  should  be  changed  to  z'  in  the  above  equation.  The  importance  and  usefulness  of  Py(z) 
comes  from  its  ability  to  reveal  some  statistical  information  about  the  underlying  process. 
Using  (2.14)  one  can  see  that 
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(2.15) 


1^1  =t«^,(n)=£W 

<fe  J,.,  ss 

and  higher  order  moments  can  be  similarly  obtained. 

Figures  2.3  and  2.4  show  a  general  and  a  detailed  version  of  the  circular  flow- 
graph  diagram  for  the  system  model  presented  in  the  previous  section.  This  gr^h  has 
a  total  of  i'+2  states.  Two  of  these  states  are  the  acquisition  state  and  the  false-alarm 
state,  where  the  former  is  absorbing  and  the  latter  is  not.  The  remaining  »=N/A  states 
are  searching  states,  where  N  is  the  PN  code  sequence  period  and  A  (taken  to  be  1/2  for 
this  system)  is  the  fraction  of  a  chip  to  be  updated  each  time  an  Hg  decision  is  made.  We 
assume  that  the  code  tracking  loop  following  the  acquisition  circuit  can  successfully  track 
the  incoming  code  phase  if  the  phase  offset  is  smaller  than  seconds.  If  we  set  7  =  1/2, 
where  yT^  is  the  timing  difference  between  the  incoming  and  the  locally  generated  PN 
waveforms  when  j=i,  then  a  total  of  four  states  will  correspond  to  ff;(alignment 
hypothesis)  as  shown  in  Figure  2.4.  The  number  of  search  states  under  Hg  is  a  v-4  or 
namely  2iV-4. 

The  search  for  the  acquisition  could  start  in  any  of  the  searching  states,  regarding 
to  the  initial  phase  of  the  incoming  PN  sequence.  Therefore,  the  initial  staitmg  state  is 
assumed  to  have  some  a  priori  distribution,  which  reflects  the  designers  confidence  about 
the  initial  relative  position  of  the  codes. 
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C.  EXPECTED  RESULTS 


Prior  to  making  a  more  detailed  analysis  this  section  deals  with  a  theoretical  best 
case  situation  to  reveal  the  relationships  between  important  parameters  for  a  serial  search 
acquisition  scheme. 

Consider  the  circular  flow-graph  diagram  of  Figure  2.3.  Let  indicate  the 
probability  that  the  system  will  reach  acquisition  while  the  codes  are  in  any  of  the  Hj 
states  and  Pj,  denote  the  probability  of  committing  an  acquisition  decision  while  the  codes 
are  in  any  Hq  states.  It  is  easy  to  see  that  P^  =  1.0  and  Pji,  =  0.0  is  the  best  case  for  this 
system.  Let  denote  the  sample  size,  i.e.,  the  number  of  chips  to  be  integrated  prior 
to  committing  a  decision.  Since  the  codes  may  take  any  initial  relative  phase,  e.g.,  as 
would  be  for  an  arbitrary  asynchronous  communication  scheme,  we  assume  a  uniform 
initial  state  prior  probability  distribution.  This  system  will  reach  the  acquisition  state  as 
soon  as  a  search  state  under  Hj  is  reached  and  wUl  never  circle  more  than  once,  so  we 
expect  to  see  the  least  possible  mean  time  to  acquire  with  the  least  variance.  It  is  easy 
to  see  that  the  only  randomness  for  the  system  is  due  to  prior  starting  state  probability 
distribution  and  the  mean  and  variance  of  the  acquisition  time  will  heavily  dq)end  on  this 
distribution.  The  mean  time  to  acquire  is 


V’TACQ  ~  ^Pi 


j'l 


Ificq 


J  2N-2 


=  TM  (N  -  -+-) 
^  ®  2  N 

“  iNT\  =  T 


(2.16) 


where  T  =  NT^  is  the  data  bit  period  and  N  is  the  PN  code  sequence  length. 
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Following  a  similar  derivation  for  the  variance  aiKl  using  M  TACQ  for 

simplicity,  we  obtain 


V  j«2  2W-2 

o«c9  •  ‘  m-mi,  -  -  Nujh 

1-1  ^  1-3 

- 

-  mf  ‘  T\ 

(2.17) 


Note  that  definite  detection  and  no  false  alarm  assumption  is  unrealistic  and  will 
not  happen  in  nature.  Nevertheless,  this  simplistic  approach  provides  some  important 
aspects  of  the  total  system. 

We  can  see  that  we  will  have  an  absolute  lower  bound  due  to  uncertainty  of  the 
phase  of  the  incoming  PN  sequence.  For  this  lower  bound  the  sample  size  constitutes  a 
key  parameter,  the  mean  acquisition  time  will  be  basicly  dependent  upon  the  sample  size 
and  the  variance  will  be  proportional  to  the  square  of  the  sample  size.  In  realistic  cases 
the  false  alarms  and  misses  will  prolong  the  average  acquisition  time  and  increase  the 
variance.  Since  the  states  under  Ho  largely  outnumbers  the  states  under  Hj,  Hq  states  will 
be  dominant  in  mean  and  variance  calculations. 

The  next  logical  question  to  ask  is  how  the  actual  detection  and  false  alarm 
probabilities  would  effect  our  performance  parameters.  The  following  sections  will 
provide  a  more  detailed  exact  analysis  for  the  overall  system  performance  and  seek 
optimal  design  parameters. 
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D.  MEAN  AND  VARIANCE  OF  THE  ACQUISITION  TIME 
1.  General 

The  flow-graph  diagram  depicted  in  Figure  2.3  provides  us  a  technique  of 
obtaining  first  and  higher  order  statistics  of  the  acquisition  time.  We  are  interested  in 
finding  the  gain  function  Pacq(z),  where 


/«i 


(2.18) 


with  Piacq(z)  denoting  the  gain  function  going  from  state  i  to  the  acquisition  state  ACQ. 
We  assume  a  uniform  a  priori  initial  phase  distribution,  i.e.,  p,  =  1/v  =  1/2N.  The 


individual  gain  functions  from  states  can  be  found  as 


PAAcq(^)  = 

PlAcq^^  ~  *  ^2^3.ac»(^)’ 

PxACqi^)  =  ^  M^ACqiZ), 

P,Acq^Z)  =  D,  -H  M  F,^,(Z). 


Note  that  these  functions  can  be  written  in  various  ways  since  we  are  dealing 
with  a  circular  graph  structure.  It  is  possible  to  write  aU  the  other  states  in  terms  of  any 
single  one  of  them.  The  form  presented  above  expresses  all  states  either  explicitly  or 
implicitly  in  terms  of  P^j,acq(^)>  the  state  which  intuitively  has  the  least  expected  time  to 
acquire  and  is  chosen  to  provide  computational  simplicity.  For  the  states  under  Hj,  M(z) 
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stands  for  the  mss  gain  function  and  stands  for  the  detectim  gain  function,  whereas 

. M^2(^)  denotes  the  process  gain  functions  for  the  states  under  Hq.  Next  we  can 

find  the  individual  gain  functions  as  follows: 


=  PJiv-\)PJip-\)z‘-"- 
=  1  - 

where  PJiv-\)  is  the  detection  probability  and  is  the  coincidence  detector 

detection  probability  for  state  v  -  i.  The  time  delay  tVi  the  sum  of  A  or  less  delay 
terms  dq)ending  on  whether  an  early  tenninatible  coincidence  d^ection  scheme  is  used 
or  not.  In  case  of  early  termination  feature  it  can  be  expressed  as 


A-B 


iU  -  E 


y*o 

«-i 

y*o 


A-j-l 

B-l 

A-B+j 

j 


(1  -p>-i)y-^-V/»'-i)y(^  7X-1 

(1  -p/p-i)r-'^*\p,ip-i)y(A  -B*i 


(2.21) 


This  equation  can  be  verified  by  examining  the  early-termination  coincidence 
detection  process  described  in  Figure  2.2.  Assuming  the  coincidence  tests  are  statistically 
independent,  the  early  termination  procedure  does  not  affect  the  coincidence  detection 
probability  which  has  a  binomial  distribution  in  both  cases.  It  can  be  written  as 
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pjy-i)  -  E  f'^ld  - 

J‘B  J 

V  ^ 


(2.22) 


Note  that  the  individual  gain  functions  for  the  other  Hj  states  can  be  found 
with  a  simple  substitution  of  the  indices.  For  the  remaining  of  the  search  states  under  Hq 
we  can  write 


M,{z)  =  NFA,{z)  +  FA^iz)  R^(z) 


(2.23) 


where 


NFA,{z)  =  (1-F^(3))z'*  +  Pj,(3)(l-P,j,(3))z'*"*', 
FA,{z)  =  P^(3)P^(3)z'‘"’  =  1  -  NFA,{z), 
R,(z)  = 


(2.24) 


In  the  above  equations  FAjfz)  and  NFAjfz)  are  the  gain  functions  for  the  false- 
alarm  and  correct  dismissal  respectively  and  R/z)  is  the  return  gain  function  for  state  3, 
which  specify  the  penalty  time  incurred  in  returning  from  the  FA  state.  The  constant  c 
is  the  number  of  bit  intervals  elapsed  before  returning  from  false  alarm  state.  Again  the 
gain  functions  for  states  4,5,6, can  be  found  with  straightforward  substitution  of 
the  indices,  whereas  the  coincidence  false  alarm  probabilities  and  time  delays  can  be 
found  from  the  equations  (2.21)  and  (2.22)  by  substituting  with  Pj^,. 

So  far  all  the  necessary  elements  constituting  the  overall  gain  function,  Pacq(^) 
are  determined.  A  simple  verification  of  the  equations  can  be  obtained  by  evaluating 
PacqO),  which  equals  to  unity  because  the  ACQ  state  is  the  only  absorbing  state  in  this 
system. 
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Following  the  stq>s  outlined  in  [Rrf.  9],  the  mean  and  the  variance  of  the 


acquisition  time  can  be  readily  obtained  by  using  the  fonnulas 


^TACQ 

„2 

°TACQ 


I 

*  Li’ 

^  d^P^iz)  ^  dP^iz) 
dz^  dz 


.  dz 


(2.25) 


2.  The  Mean  Acquisition  Time 

Now  we  will  carry  out  some  algd>ra  starting  from  the  overall  gain  function 
Pacq(z)  of  equation  (2.19).  The  following  symmetry  conditions  always  hold  if  complete 
random  sequences  were  used.  They  are  also  valid  for  m-sequences  with  long  periods. 


KM)  A  M  (z)  =  M,(z) 

MJz)  A  Af,.,(z)  =  Af,(z) 

A  D,(z)  =  D,(z)  (2.2«) 

D^(z)  A  D  .,(z)  =  Dj(z) 

A  M,iz)  =  M,iz)  =...=  M^.,iz) 

Note  also,  with  ACQ  being  the  only  absorbing  state  =  1 .0,  i  =  1 , . . .  ,1'  and 

A^(l)  =  1 .0.  Using  these  results  and  carrying  out  the  algebra  we  can  write 


pLM)  =  (*'-4)<(l)  +  PU^il), 

P^qil)  =  D^l)  +  a4(1)  +  Af,2(l)P;^(l). 
pUH)  =  D'liD  *  +  M,,il)P;^il), 

nU(l)  =  D'jH)  *  Mj,(l)  ^  M,,{l)PU,il) 

where 
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<(1)  = 

Dkl)  = 

and  and  follows  with  an  indices  substitution. 

The  other  states  under  Hq  can  be  similarly  found  as  below, 

Pia,(Xi  =  (v-i-l)A#ji(l)  +  ,  i  =  4.5 . v-2 


(2.28) 


(2.29) 


where 


^  ti.P^  *  cNT^Pj,P^. 


(2.30) 


Finally  for  the  state  v-l  we  can  write 


^jzCl)  ,  ,  , 

+  - — — - {2MJMJ>^+2MJmJD., 

+  - 

(1-M  ^(1))^  dJ  d2  J2  dl  d2  dl 

*  "  ■"  ,231) 


After  r^lacing  v  by  2N,  the  expression  for  the  mean  acquisition  time  for  our 
system  can  be  written  in  the  following  form 
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(2.32) 


l*TACa  ~ 

/•I 

=  JL[(2N*-7iV+6)Af;(l)+(2iV-3)P;.,^(l)+P/^(l) 

2N 

3.  Hie  Variance  of  Hie  Acquisition  Hme 

Equation  (2.25)  dqiicts  that  an  expression  for  the  variance  can  be  obtained  by 
taking  the  second  derivative  of  P^(z).  Using  the  same  symmetiy  conditions  shown  in 
equation  (2.26)  we  can  write 


PiU^)  =  ((i'-4)i<  +  (r-4)(r-5)(w;)^ 
+  (2i.-8)Af>;-,^  +  Piix^)il), 
P2U^)  =  (PH2  +  MH2  *  23/^3^  + 
P'xU^)  =  (P'Jt  *  M41  +  2Af>3'^  + 
PlUh  =  (P'J,  *  <  *  23f>U  ^ 


(2.33) 


and  for  the  other  states  under  Hq 


p!U^)  *  ^  (r-l-l)(r-/-2X3/;)* 

+  2(r-i-l)3/>;.,^  +  P^'.^^{\),  i  =  4,5 . r-2. 


(2.34) 


Hie  gain  function  P^i^(z)  has  the  form  A(z)fB(z)  and  can  be  differentiated 


using  normal  techniques. 


(2.35) 


where 
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A(z)  =  D^iz)  *  MJzyD,jiz)  *  MJ^z)M,jiz)D,,iz)*  M„iz)Mjiz)DJz), 
B(z)  =  1  -  MjizWji^zW/^iz). 


(2.36) 


The  rest  follows  with  a  tedious  but  simple  bookkeq)ing  exercise.  The  resulting 
equations  which  will  be  used  to  calculate  the  variance  of  the  acquisition  time  are  rather 
messy  and  we  do  not  show  them  here,  but  they  can  be  easily  put  in  a  computer  program 
and  an  acquisition  time  analysis  may  follow. 
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m.  ACQUISITION  TIME  ANALYSIS 


A.  INTRODUCTION 

The  decision  processor  of  the  DS-SS  receiver  dq>icted  in  Figure  2.1  has  the  task 
of  resolving  between  the  hypotheses  given  by  equation  (2.7).  This  task  will  be 
accomplished  using  the  test  statistics  Y„.  The  decision  processor^  can  be  implemented  in 
various  schemes  based  on  the  likelihood  ratio  of  equation  (2.8). 

In  this  chapter,  we  will  study  three  decision  schemes;  a  fixed  sample  size  (FSS) 
test  scheme,  a  sequential  probability  ratio  test  (SPRT)  scheme  and  a  truncated  SPRT 
scheme.  These  schemes  are  previously  studied  in  [Ref.  7]  and  [Ref.  10]  for  a  single 
search  state  cell  basis  to  obtain  the  average  test  Imigths  for  SPRT  and  TSPRT  and 
comparing  them  to  that  of  the  FSS  test.  The  flow-grcq>h  technique  presented  in  previous 
chapter  is  applicable  to  the  all  decision  schemes  using  a  serial  search  algorithm,  hence 
can  be  adequately  used  to  study  some  important  design  considerations  in  terms  of  the 
performance  parameters,  namely  the  mean  and  the  yariance  of  the  acquisition  time. 

B.  FIXED  SAMPLE  SIZE  SCHEME 

In  a  fixed  sample  size  (FSS)  decision  scheme,  the  length  of  integration  is  fixed  and 
a  decision  is  made  based  on  the  resulting  test  statistics,  Y„.  If  the  integration  length  is 
taken  as  from  t  =  Oto  t  =  MT„  the  test  can  be  described  by 

The  decision  processor  design  methods  discussed  in  this  chapter  are  adapted  from  [Ref.  7]. 
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(3.1) 


f  ^  T  say 
I  <  T  say 


where  t  is  the  threshold.  Since  the  likelihood  ratio  of  equation  (2.8)  is  a  monotonically 
increasing  function  of  the  variable  equivalently  the  test  can  be  written  as 


^  t'  =  A«(t)  say  /f, 
<  t'  =  Am(t)  say 


where  A„  '  (.)  is  the  inverse  function  of  AJ,.)  of  equation  (2.8). 

The  cumulative  distribution  function  (cdf)  of  Y„  can  be  written  in  terms  of  the 
Marcum  Q  function  as 


Fk.W  =  - 1  - «  y.  ^  0 


where  the  Q  function  is  defined  as  [Ref.  11] 


<?(C,5)  =  /*  dx. 


An  iterative  algorithm  for  calculating  the  Q  function  can  be  found  in  [Ref.  12]. 
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1.  De^gn  of  Decision  Panuneters 

Using  equation  (3.3)  we  can  write  the  false  alarm  and  the  miss  probabilities 
for  the  FSS  test  as 


miss 


(K 

^iV.O 
— T"  » 

4) 

Ou 

Ou 

1  - 

(K 

2 

> 

1 

Ou 

(3.5) 


4)- 

<^A# 


From  these  equations,  the  values  of  r'  and  M  can  be  obtained  by  iteratively  solving 
the  equations  simultaneously,  such  that  Pj^  ^  a  and  ^mlss  ^  ^  -  P,  where  a  and  jS  are 
the  desired  false  alarm  and  detection  probabilities  req)ectively. 

The  nominal  worst  case  values  of  (KJo^)  are  derived  in  Chapter  n  by  equation 
(2.13)  and  repeated  below. 
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(3.6) 


Considering  the  structure  of  flow  gr^h  dqjicted  in  Figure  2.3  of  Chapter  II  we 
will  need  the  ratio  for  the  case  of  |  j+y  |  =1.5.  Recall  from  equation  (2.7) 

that  for  I  j +7  I  ^  .5  we  have  the  hypothesis  and  for  |  j+y  |  ^  2.0  we  have  the 
hypothesis  Hq.  Therefore  (\,/20  ratio  for  |  j+7  |  =1.5,  denoted  by  ()^  ^rJ2a^  can 
be  £q}proximated  by 
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(3.7) 


^it,>/2  _  1  \.l  ^  2  \ft 

2al  3  2al  3  2al 
using  a  first  order  lii^ar  i^roximation. 

2.  Mean  and  Variance  Analysis  of  The  FSS  Schone 

We  are  now  ready  to  use  the  expressions  derived  for  the  mean  and  variance 
of  the  acquisition  time.  Note  that  the  required  parameters  are  Pji  ^  /3,  =  a,  all  the 

time  delays  for  initial  tests  are  MT^  and  the  time  delays  for  coincidence  tests  can  be 
calculated  using  equation  (2.21).  Finally, 


a.  Optimal  Choices  of  The  Design  Parameters 

We  have  seen  that  specifying  the  desired  false  alarm  and  detection 
probabilities  suffices  to  design  a  FSS  decision  scheme.  In  previous  studies  these 
parameters  are  forced  to  be  arbitrarily  chosen  obeying  the  intuition  that  a  should  be 
reasonably  close  to  zero  while  should  be  reasonably  close  to  unity.  Knowledge  of  the 
mean  and  the  variance  of  the  resulting  acquisition  time  provides  us  a  means  of  optimal 
selection  of  these  parameters  with  respect  to  these  performance  parameters. 

Unfortunately,  equation  (3.5)  dqpicts  that  a  and  &  are  related  with  the 
design  parameters,  the  sample  size  and  the  threshold  in  terms  of  Marcum  Q  function. 
The  values  of  interest  does  not  allow  us  to  make  the  well  known  exponential 
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approximations  of  the  Marcum  Q  function,  hence  we  have  to  depend  on  the  numerical 
results  rather  than  a  more  desirable  rigorous  analytical  analysis. 

Our  argument  will  now  continue  with  the  illustration  of  some  examples. 
As  the  fu^t  example  we  employ  a  FSS,  PN  code  serial  search  acquisition  scheme  with 
an  early  terminatable  coincidence  detector.  Using  the  notation  of  Chapter  n  we  assume 
7  =  .5,  A  =  .5  for  this  system.  The  coincidence  detector  parameters  are  A  =  4  and  B 
=  2.  The  design  chip  SNR  =  -10  dB  and  the  amount  of  time  elapsed  for  returning  from 
the  FA  states  (penalty  time)  is  set  to  be  S  bit  periods. 

Figures  3.1  and  3.2  show  the  graphs  for  normalized  Mtacq  and  o^tacq 
versus  various  values  of  desired  false  alarm  (a)  and  detection  (j3)  probabilities  for  this 
system  respectively.  The  plots  have  the  following  common  properties: 

•  The  rate  of  the  change  of  the  values  of  with  a  are  much  greater  than  the  variations 
due  to  /3. 

•  Both  curves  has  a  minimum  point  for  a  certain  a  and  combinations. 

The  minimum  ftTAcq  is  111.46  T and  occurs  at  a =.06,  /8=.69,  whereas 
the  minimum  o^tacq  is  6.0e3  T  ^  and  occurs  at  a  =  .06,  j8  =.8. 

Let  us  consider  the  case  of  minimum  Mtacq  carry  out  a  simple 
comparison  with  the  theoretical  values  computed  in  Section  C  of  Chapter  H.  The  fixed 
sample  size,  Mf„  calculated  by  equation  (3.5)  for  these  values  of  a  and  |8  is  70,  and 
recall  from  equation  (2.17)  that  70  bit  periods  is  the  lowest  mean  time  to  acquire  if  it 
were  possible  to  achieve  definite  detection  and  no  false  alarm  hypothesis  with  this  sample 
size  value.  On  the  other  hand,  the  real  values  are  0.6967,  Pj2=  0.2846,  Pji,=  0.06, 


the  coincidence  test  detection  and  false  alaim  probabilities  from  equation  (2.22)  are  = 
0.9138,  Pcj2—  0.3213  and  P^=  0.0199.  The  coincidence  detector  time  delays  from 
equation  (2.21)  are  /j7'=189.5,  /j2'=234.9  and  tj^'=22Q.%  chip  periods.  The  effective 
sample  size  per  state,  i.e.,  M^'(l),  for  Hq  states  can  be  calculated  using  equation  (2.30) 
to  be  89.36.  Therefore,  considering  a  /^=0.06,  the  system  would  have  reached 
acquisition  on  the  average  89.36  bit  periods  if  it  were  quaranteed  to  acquire  as  soon  as 
it  reaches  a  state  under  //,.  However,  the  actual  mean  acquisition  time  came  out  to  be 
111.46  bit  periods.  Hi  states  and  additional  sweeps  of  the  uncertainty  region  prolonged 
the  mean  acquisition  time  to  1 1 1 .46,  only  22. 1  bit  periods.  As  we  remarked  earlier  the 
states  under  Hq  dominates  the  resulting  mean  and  variance  of  the  acquisition  time.  This 
simple  example  and  the  greater  rate  of  change  with  respect  to  a  verifies  this  hypothesis. 

Surprisingly,  we  observe  that  the  desired  detection  probability,  /8  which 
minimizes  either  the  mean  or  the  variance  of  the  acquisition  time,  takes  moderate  values, 
opposing  the  intuitive  thought  that  it  should  be  as  close  to  unity  as  possible.  We  will 
explain  this  phenomenon  by  stressing  the  importance  of  the  sample  size  per  decision.  We 
have  seen  in  Section  C  of  Chapter  n  that  the  average  acquisition  time  heavily  dq)end  on 
the  sample  size.  Therefore,  it  should  be  kept  as  low  as  possible  to  minimize  the 
acquisition  time.  For  a  fixed  dwell  acquisition  scheme  such  as  FSS,  the  sample  size  per 
decision  is  fixed  a  priori  and  we  denote  it  The  plot  of  the  Mp,  versus  a  and  /3  is 
given  in  Figure  3.3. 
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Figure  3.1  Average  acquisition  time  for  various  a  and  j8  combinations. 


Figure  3.2  Variance  of  the  acquisition  time  for  various  a  and  fi  combinations. 
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Figure  3.3  The  fixed  sample  size  (Mfss)  versus  various  a  and  combinations. 

This  graph  has  an  interesting  property  that  along  the  main  diagonal  the 
values  are  "almost”  symmetrical,  that  is,  for  a  fixed  a,  Afjs,  increases  exponentially  with 
increasing  /3,  and  for  a  fixed  decreasing  a  has  "almost"  the  same  effect.  In  fact  the  rate 
of  change  with  respect  to  is  slightly  greater.  Choosing  a  small  a  and  large  /3,  requires 
a  large  which  in  turn  results  in  a  large  mean  and  variance  of  the  acquisition  time. 
Otherwise,  false  alarms  and  additional  sweqjs  of  the  uncertainty  region  increase  the  same 
parameters.  Therefore,  an  optimum  acquisition  time  can  be  achieved  in  terms  of  the 
design  parameters  a  and  /3.  Since  the  states  under  bear  far  more  weight  than  the  states 
under  ffj  and  false  alarms  are  costly,  minimal  a  turns  out  to  be  a  small  value  (close  to 
zero),  while  selecting  a  large  increases  the  Afj^,  excessively  with  respect  to  its 
contribution,  minimal  jS  takes  moderate  values. 
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The  second  example  is  chosen  to  illustrate  the  potential  benefits  of  using 
a  coincidence  detector  in  the  acquisition  system.  The  only  difference  is  that  we  have  s^ 
the  values  of  A  and  B  both  to  be  zero,  and  therefore  no  coincidence  detector  is  used.  The 
minimum  value  of  the  average  acquisition  time  is  now  increased  to  171.9877  T  for  the 
values  of  a =0.0046,  /3=0.63.  The  minimum  value  of  the  variance  is  also  increased  to 
13702  for  a  =  0.0051,  /3  =  0.79. 

Our  ailment  for  the  design  parameters  completely  fit  in  this  system  as 
well.  Note  that  minimal  a  turns  out  to  be  ly^roximately  10  times  smaller  than  the 
previous  example,  because  coincidence  detector  by  design  decreases  the  effective  false 
alarm  probability.  We  will  investigate  this  matter  further  in  detail. 

Our  discussion  so  far  targeted  the  two  performance  parameters,  i.e.,  the 
mean  and  the  variance  of  the  acquisition  time  sq)arately.  We  have  seen  from  the 
examples  that  both  parameters  are  minimized  by  close  but  distinct  choices  of  design 
parameters  a  and  jS  and  yet  no  suggestion  has  been  made  which  pair  should  be  used  in 
the  design.  Before  further  discussion,  we  summarize  our  results  for  the  design  examples 
in  Table  3.1. 

The  Tchebycheff  Inequality  [Ref.  13]  is  a  well  known  theorem  of 
probability  theory  and  is  a  measure  of  the  concentration  of  an  RV  near  its  mean  and 
its  variance  o^.  It  can  be  formulated  as 

^  —  ,V  e  >  0.  (3.9) 
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TABLE  3.1  The  system  values  for  the  design  examples  of  FSS  test. 


Minimum  Mtacq  Parameters 

Minimum  o^tacq  Parameters  | 

M) 

Mtacq 

U^TACQ 

(a./3) 

Mtacq 

TACQ 

System 

with 

ETCD 

(0.06, 

0.69) 

111.46 

7014.0 

(0.06, 

0.80) 

119.08 

6002.0 

System 

with 

No  CD 

171.99 

16572.0 

(0.0051, 

0.79) 

183.63 

13702.0 

The  proof  of  this  theorem  is  simple  and  can  be  found  in  [Ref.  13].  Note 
that  even  the  bound  in  (3.9)  is  shown  to  be  not  so  tight  for  specific  densities,  i.e.,  normal 
densities,  we  could  still  use  it  conceptually  to  find  a  criteria  for  choosing  optimal  design 
parameters. 

We  remark  that  the  probability  density  function  of  the  acquisition  time 
be  positive  definite,  that  is  "^ACQ  ^  0'  We  are  interested  in  Hnding  an  upper  bound  for 
the  probability  for  the  cases  "^ACQ  is  greater  than  /itacq  up  to  several  order  of  magnitude. 
Let  denote  the  maximum  "acceptable"  acquisition  time  for  a  system  of  this  nature. 
Setting  e  =  T^-  (trACQ  in  (3-9),  we  obtain 


<  P{\ 

^ACQ  ~  I^TACqI 


^  ^VBM3i~PTAC<^  ^ 


2 

<^TACQ 


PtaCc)^ 


(3.10) 


To  clarify  a  possible  ambiguity,  the  reader  should  note  that  we  are  not 
actually  interested  in  determining  the  probability  that  the  system  will  reach  acquisition 
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in  less  than  seconds,  although  equation  (3.10)  dq>icts  a  pessimistic  and  quite 
possibly  a  loose  bound  for  this  probability.  The  usefulness  of  the  Thebycheff’s  Inequality 
comes  from  the  fact  that  it  can  be  used  for  any  random  variable  if  the  first  and  second 
order  statistics  are  known.  Therefore,  it  can  be  adequately  used  for  our  case,  to  make 
an  even  comparison  between  the  multiple  choices  with  which  we  are  confronted  and  to 
help  us  to  narrow  our  attention  to  a  unified  direction. 

The  right-most  expression  in  (3.10)  is  the  bound  for  the  upper  tail 
probability  which  we  want  to  minimize.  To  be  consistent  with  the  notation  we  have  been 
using,  we  will  denote  the  parameters  that  minimizes  the  right-most  expression  using  an 
asterisk  (*)  following  the  relevant  parameter  to  demonstrate  that  the  parameter  is 
"optimal"  in  this  context.  The  numerical  results  for  various  T^  suggests  the  a*  and  /S* 
combinations  lay  between  the  cases  of  minimum  Mtacq  minimum  o^tacq-  I**  f^^ct  it 

is  easy  to  see  from  equation  (3. 10)  that  for  the  case  of  >  Mtacq  minimum  o^tacq 

choices  minimize  the  right-most  expression  in  (3.10).  Figure  3.4  dq)icts  this  behavior 
of  a*  and  /3*  graphically  for  the  design  examples  of  this  section  from  Table  3.1.  System 
1  corresponds  to  the  design  example  which  uses  an  early  terminatable  coincidence 
detector,  whereas  system  2  is  the  example  without  the  coincidence  detector  in  this  figure. 
The  X  -  axis  assumes  various  values  of  T„,g/r  normalized  by  the  sequence  length  N. 
Note  that  a*  values  remain  essentially  unchanged  for  the  entire  range,  which  verifies  our 
earlier  observation  about  the  dominance  of  the  Hg  states.  Figure  3.5  shows  the 
corresponding  values  of  the  normalized  /**tacq  <^*tacq  of  the  same  examples  for  the 
same  range  of  /NT.  Note  that  as  increases  m’tacq  values  increase  and  o^Vacq 
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Figure  3.4  Optimal  a  and  /3  values  with  respect  to  the  Thychebyceffs  Inequalit} 


Figure  3.5  Optimal  performance  parameters  with  respect  to  the  Thychebycheffs 
Inequality. 
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values  decrease,  but  the  change  in  either  of  these  is  not  so  drastical.  This  is  because  these 
parameters  are  essentially  determined  by  a*,  which  remains  constant  for  the  whole  range. 

Had  the  value  of  been  specified  by  the  system  designer,  we  could 
use  that  value  in  our  equations  for  further  discussion.  Since  we  (k)  not  have  a 
predetermined  value,  we  will  narrow  our  focus  to  the  minimum  o^tacq  choices  for  the 
remainder  of  this  thesis.  We  will  also  present  the  results  for  minimum  ftrAcg  choices 
where  appropriate. 

b.  The  Eiffeets  of  The  Penalty  Time 

The  penalty  time  for  our  system  was  defined  as  the  amount  of  bit 
intervals  elapsed  for  returning  from  the  FA  state.  It  is  essentially  determined  by  the 
external  tracking  circuit  as  a  measure  of  the  quality  of  the  scheme.  Recall  from  Ch2q>ter 
n  that  the  penalty  time  was  expressed  as  cNT^,  where  c  denoted  the  number  of  data  bits, 
N  was  the  sequence  length  and  T^  was  the  chip  period.  Using  the  flow  graph  technique 
and  selecting  the  minimum  variance  choices  we  have  calculated  a*,  |8*,  and  normalized 
m’tacq  and  o^’tacxj  while  the  c  parameter  is  varied  from  1  bit  period  to  1000  bit  periods. 
Figure  3.6  depicts  the  results  for  the  system  of  example  1  from  Table  3.1. 

Increased  penalty  time  prolonged  the  average  acquisition  time  and 
increased  the  variance  as  naturally  expected.  As  far  as  the  design  parameters  are 
concerned,  we  see  that  /3*  is  not  affected,  since  the  penalty  only  occurs  at  Hq  states.  The 
a*  is  not  affected  very  much  either,  if  the  variation  of  the  penalty  time  remains  in  close 
vicinity  of  its  design  value.  The  system,  however,  requires  a  much  lower  design  a  to 
compensate,  for  an  excessively  large  penalty  time. 
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e  (No  of  Peftolty  BHs) 

Figure  3.6  The  effects  of  penalty  time  on  the  design  parameters. 

Using  the  design  values  shown  in  Figure  3.6,  we  can  now  examine  the 
system  behavior  for  the  cases  where  actual  penalty  time  differs  than  specified  in  the 
design.  We  have  considered  three  design  values  of  the  penalty  time,  where  c  =  S,  SO  and 
SOO  bit  periods.  The  system  is  constructed  using  the  corresponding  a  and  for  each  case 
from  Figure  (3.6).  The  actual  penalty  time  is  again  varied  from  1  to  1000  bit  periods. 
The  variation  of  o^tacq  changing  penalty  time  is  given  in  Figure  3.7.  We  saw 

that  if  the  system  spent  more  time  than  specified  in  the  design  to  recover  from  a  false 
alarm  state,  this  would  degrade  the  system  performance.  This  degradation  would  effect 
the  system  more  severely,  especially  when  the  specified  penalty  time  is  small  and  actual 
penalty  takes  considerably  longer  time.  Consider  our  example,  the  system  designed  using 
c  =  S,  results  in  a  variance  value  1.95  times  higher  than  the  system  designed  using  c  = 
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so  when  the  actual  pmalty  time  (Ca)  is  SO  bit  periods,  whereas  the  seomd  syston  has 
only  1 .2  times  greater  variance  value  than  the  first  system  at  c^  =  S  and  the  degradation 
gets  worse  as  Ca  increases.  Although,  we  have  shown  here  some  extreme  cases,  this 
behavior  of  the  system  is  still  valid  for  a  smaller  scale  of  the  change  in  the  penalty  time. 
Therefore,  we  can  say  that  if  some  uncertainty  exists  about  the  penalty  time,  it  is  safer 
to  design  the  system  using  highest  possible  c  value. 


ACTU«.  PENALTY  TIME 

Figure  3.7  The  variation  of  o^tacq  with  the  actual  penalty  time. 

c.  The  Channel  Mismatch  Problem 

We  will  now  examine  a  more  serious  and  practically  impoitant  design 
problem  of  the  acquisition  system.  Recall  from  equation  (2.11)  that  the  chip  signal  to 
noise  ratio  (SNR)  of  the  incoming  signal  had  to  be  specified  prior  to  designing  the 
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system.  In  practice,  for  the  receiver’s  site  there  is  no  way  of  knowing  the  exact  SNR 
other  than  obtaining  this  value  by  field  measurements  and  this  effective  SNR  value  may 
be  different  than  ^lecifications  for  a  vari^  of  reasons,  i.e.,  the  receiver’s  location, 
channel  medium,  fading. 

Using  equation  (2.11)  we  can  show  that  the  bit  signal  to  noise  ratio, 
denoted  by  SNR^  can  be  expressed  as  SNRi,=  (Aq^  T  /  2N^  =  N  SNR.  Here  N  is  the 
sequence  length  and  T  is  the  bit  period.  So  far,  for  our  examples  we  fixed  the  per-chip 
SNR  as  -10  dB,  which  corresponds  to  an  SNRf,  of  20  dB.  Table  3.2  dq>icts  various  design 
SNR  values  and  the  resulting  /i*TACQ  ^<1  <’^*tacq  using  minimum  variance  choices  for 
the  system  of  example  1. 

We  will  define  the  term  "channel  mismatch"  as  the  response  of  the 
receiver  to  the  incoming  signals  that  have  different  SNR  than  specified  in  the  initial 
design.  Naturally,  under  occurrence  of  such  an  event  some  performance  degradation 
should  be  expected.  This  degradation  wiU  luq>pen  for  the  stronger  signals,  as  well  as  for 
the  weaker  signals,  because  lower  design  SNR  will  yield  a  lower  threshold  which  can 
result  an  increase  in  the  number  of  false  alarms  for  stronger  signals. 

Using  the  flow  graph  technique  we  present  the  variations  of  resulting 
mean  and  the  variance  of  the  acquisition  time  with  respect  to  the  effective  SNR  of  the 
incoming  signal.  The  system  of  example  1  is  used  for  three  sqiarate  design  SNR  values 
(-13,  -10  and  -3  dB)  from  Table  3.2.  The  effective  SNR  is  changed  from  -20  dB  to  0  dB 
for  all  three  cases  and  the  plots  for  Mtacq  and  o^tacq  are  given  in  Figures  3.8  and  3.9 
respectively. 


40 


TABLE  3.2  The  system  values  for  various  design  SNR  for  the  acquisition  system. 


SNR(dB) 

SNRb(dB) 

• 

a 

/S’ 

• 

M  TACQ 

^  TACQ 

-20.0 

10.0 

0.12 

0.75 

927.0272 

4.2269e5 

-13.0 

17.0 

0.07 

0.79 

219.1580 

2.0658e4 

-10.0 

20.0 

0.06 

0.80 

119.0972 

6.0020e3 

-  3.0 

27.0 

0.04 

0.81 

39.9832 

4.4971e2 

0.0 

30.0 

0.03 

0.81 

20.4277 

1.6984e2 

For  all  of  the  three  cases,  we  can  say  that  the  degradation  caused  by  the 
stronger  signals  are  much  less  than  by  weaker  signals.  For  example;  for  the  case  of 
design  SNR  =  -10  dB,  a  10  dB  increase  in  the  effective  SNR  increases  /ixAco  by  a  factor 
of  3,  whereas  a  -10  dB  decrease  increase  the  same  parameter  by  a  factor  of  100.  Also, 
one  can  observe  that  in  the  close  vicinity  of  the  design  SNR  slightly  stronger  signals  yield 
a  better  performance. 

Therefore,  if  the  system  designer  has  some  uncertainty  about  what  the 
received  SNR  would  be,  it  is  better  to  design  the  system  using  the  lowest  possible  SNR 
to  counterbalance  the  channel  mismatch  problem. 

3.  The  Coincidence  Detector  Parameters 

In  Chapter  n  we  have  proposed  a  majority  logic  verification  scheme  for  our 
system  to  counter  balance  the  effects  of  costly  false  alarms.  Briefly,  the  system  will 
enforce  A  or  less  additional  tests  in  case  of  an  Hj  decision  to  verify  the  correctness.  If 
B  of  the  additional  tests  are  in  favor  of  the  original  decision  the  system  decides  Hj.  We 
also  assumed  that  a  counter  is  available  and  the  test  can  be  early  terminatable  once  B 
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10-*  10-'  10« 
EFFECTIVE  PER-CHIP  SNR 

Figure  3.8  The  variation  of  /xtacq  with  the  effective  per-chip  SNR. 


Figure  3.9  The  variation  of  o^tacq  with  the  effective  per-chip  SNR. 
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favorable  tests  are  accumulated  or  there  is  no  chance  that  B  favorable  tests  can  be 
obtainable.  In  this  section  we  will  investigate  what  good  choices  of  the  param^rs  A  and 
B  should  be,  with  the  help  of  the  flow  graph  technique. 

We  have  calculated  minimum  o^tacq  and  corresponding  a  and  jS*  values  for 
several  coincidence  detector  (CD)  parameter  pairs  (A,B)  for  the  system  of  example  1 , 
both  for  with  and  without  early  termination  feature.  We  present  the  results  in  Figure  3. 10 
graphically. 

The  best  case  for  the  early  terminatable  coincidence  detector  (ETCD)  occurs 
for  the  (A,B)  pair  (5,3)  resulting  a  minimum  variance  of  5715.8  7^  with  a*  =  .09  and 
/3*  =  .82  and  the  corresponding  /xtacq  =  117.07  T.  This  minimum  variance  value  is  the 
overall  best  case  we  have  achieved  for  the  system  with  SNR=^-10  dB  so  far.  For  the  non- 
early  terminatable  coincidence  detector  (NETCD)  we  obtain  the  minimum  o^tacq  = 
6486.4  at  (4,2)  with  a*  =  0.056  and  /3*  =  0.8,  and  the  corresponding  123.94  T. 

We  can  remark  that  choices  of  optimal  design  parameters  does  not  change 
very  much  regarding  to  the  coincidence  detection  scheme.  One  can  also  see  that  for  this 
minimum  variance  cases  the  improvement  achieved  using  an  early  terminatable 
coincidence  detector  is  approximately  5  %  in  the  mean  12  %  in  the  variance. 

Although  we  show  only  a  portion  of  the  results,  we  have  also  studied  the 
cases  of  7  ^^4  ^9  and  the  results  are  omitted,  because  none  beats  the  minimum  choices 
presented  above.  In  fact,  it  is  logical  to  conclude  that  choosing  A  too  large  is 
unacceptable  for  both  schemes,  although  using  a  ETCD  increases  the  range  a  little. 
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Figure  3.10  Minimum  o^tacq  versus  various  coincidence  detector 

parameters. 

From  Figure  3.10,  we  can  conjecture  that  once  A  is  fixed  the  optimal  choice 
of  B  is  around  \_A/2\  .  We  can  qualitatively  verify  this  hypothesis  by  the  following 
argument.  Recall  from  Chapter  n  that  the  coincidence  detection  or  false  alarm  probability 
is  given  by  the  cumulative  binomial  distribution  of  the  form 


=  p(B^  X  ^  f^l  pK\-pY-^ 

j-B 


(3.11) 


where  the  random  variable  X  represents  the  event  that  an  acquisition  decision  (either 
genuine  or  falsely  )  is  made  and  p  denotes  the  detection  or  false  alarm  probability  for  a 
single  test.  This  distribution  has  the  mean  pA  and  the  variance  p(l-p)A.  The  value  of  P^ 
equals  to  unity  for  B=0  regardless  of  the  value  of  p.  Considering  (3. 1 1)  one  can  see  that 
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for  a  fixed  A,  a  small  p,  i.e.,  a  typical  false  alarm  probability,  results  in  a  small  mean 
and  variance,  therefore,  decreases  rapidly  with  the  increasing  B.  On  the  other  hand 
a  large  p,  i.e. ,  a  typical  detection  probability,  results  a  large  mean  (close  to  A)  and 
variance,  therefore,  rapid  deteoration  in  occurs  for  large  values  of  B.  The  system 
does  not  know  in  advance  whether  an  initial  Hj  decision  will  result  in  FA  otACQ  state, 
therefore,  uses  the  same  CD  parameters  for  both  cases.  We  want  a  large  B  to  pull  down 
the  P^  as  low  as  possible,  while  a  small  B  will  make  as  high  as  possible.  Therefore 
an  optimum  B  exists  in  the  interval  [0  A].  Considering  all  the  other  factors  effecting  the 
resulting  /xxacq  ^ind  o^tacq.  it  happens  to  be  around  \A/2\  . 

4.  Simulation  Results 

Our  study  on  the  FSS  serial  search  acquisition  scheme  was  based  on 
approximations  and  analytical  formulations.  In  this  section  we  present  simulation  results 
for  some  particular  cases  to  verify  the  validity  of  our  assumptions  and  analyses. 

We  have  used  an  m-sequence  of  length  1023  with  the  primitive  generator 
polynomial  1  +  x^  +  x*  +  x*  +  x'°.  The  value  of  y  used  for  the  design  is  0.5.  The 
update  fraction,  A,  of  the  local  m-sequence  is  set  to  be  1/2,  half  the  chip  period.  The 
initial  state  a  priori  distribution  is  assumed  to  be  uniform.  The  early  terminatable 
coincidence  detector  parameters  are  taken  as  A=  A  and  B=2. 

The  simulations  are  performed  using  the  three  different  per-chip  design  SNR 
values  from  Table  3.2  (-13dB,-10dB  and  -3dB)  for  10(X)  runs.  The  sample  and  theoretical 
values  of  ^tacq  and  o^tacq  are  given  in  Table  3.3. 
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TABLE  3.3  The  simulation  results  for  various  per-chip  SNR  . 


1  SNR  (dB) 

Sample 

Mtacq 

Theoretical 

Mtacq 

Sample 

O^TACQ 

Theoretical 

<^TACQ 

-13.0 

217.10 

219.16 

19212.0 

20658.0 

-10.0 

116.33 

119.09 

5494.6 

6002.0 

-3.0 

36.82 

39.98 

574.7 

449.7 

Finally,  we  present  the  histogram  for  the  case  of  per-chip  SNR  =  -10  dB  to 
give  an  idea  about  the  shape  of  the  probability  distribution  function  fTacq(tacg)  in  Figure 
3.11. 
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Figure  3.11  Histogram  for  the  simulation  of  the  FSS  scheme 
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C.  SEQUENTIAL  AND  TRUNCATED  SEQUENTIAL  PROBABILITY  RATIO 
TESTS 

The  sequential  probability  ratio  test  (SPRT)  for  the  acquisition  scheme  of  Chapter 
n,  is  obtained  by  testing  the  likelihood  ratio  of  (2.8)  against  two  thresholds  for 
n = 1 ,2,3. . .  until  one  of  the  thresholds  is  exceeded.  The  length  of  the  integration  increases 
by  one  chip  each  time  n  increases  by  one.  The  resolvement  between  the  hypotheses  Hq 
and  Hi  is  achieved  by  checking  which  threshold  value  is  reached  fu^t.  The  test  can  be 
written  as 


SPRT: 


A„(y„) 


^  A  say 

^  B  say  Hq 

otherwise,  continue  to  next  n 


(3.12) 


Note  that  Y„  and  Y„+]  are  not  independent  and  the  likelihood  ratio  AJyJ  can  not 
be  written  as  the  sum  of  independent  random  variables.  However,  Type  I  error  (false 
alarm  probability)  and  Type  n  error  (miss  probability)  can  still  be  related  through  Wald’s 
inequalities[Ref.  14], 


and  (3.13) 

a  1  -a 

where  a  and  (= 1-jS)  are  the  resulting  false-alarm  and  miss  probabilities,  respectively. 
If  the  excess  over  the  boundary  is  small  when  the  test  terminates,  which  is  the  case  when 
the  average  test  length  is  large,  the  inequalities  in  (3.13)  can  be  approximated  by 
equalities.  Since  l-/3„<  1  and  l-a<  1,  we  can  write 
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and 


/3.  <  B 


(3.14) 


a  < 


which  is  useful  when  strict  bounds  on  the  errors  are  desired.  These  bounds  are 
tight  when  the  errors  and  the  signal-to-noise  ratio  are  small. 

Using  the  monotonicity  of  the  likelihood  ratio,  the  test  in  (3.12)  can  be  written  in 
a  more  practical  form  for  real  time  implementations,  as 


k  Ain)  =  A«(/f)  say 
S  Bin)  =  a;'(B)  say 

otherwise,  continue  to  next  n 

where  the  thresholds  can  be  precomputed  as  functions  of  n. 

The  SPRT  has  been  shown  to  be  optimum,  in  the  sense  of  minimizing  the  average 
test  length  under  and  Hj,  only  when  the  samples  are  i.i.d.[Ref.  6].  But  for  our 
situation  the  samples  are  non-i.i.d.  and  no  optimality  has  been  established.  However,  we 
expect  that  the  performance  is  still  good  and  much  better  than  the  FSS  test. 

The  SPRT  test  has  a  drawback  that  for  certain  PN  sequence  phase  disparities,  it 
may  take  an  excessively  long  time  to  resolve  between  the  two  hypotheses.  Such  an 
incident  is  especially  likely  when  the  initial  phase  difference  between  the  two  PN 
sequences  fall  in  the  state  between  Hq  and  Hj  as  shown  in  (2.7).  To  avoid  a  lengthy  test, 
a  truncated  SPRT  (TSPRT)  test  was  proposed  in  [Ref  7],  which  imposed  an  upper  bound 


on  the  test  length.  The  test  is  described  by 
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TSPRT: 


if  nt^A, 


< 


^  A  say 

^  B  say 

otherwise,  continue 


(3.1Q 


K<y^ 


if  n=A, 


A„(yJ  =  A,(y^ 


s  f  say 

<  T  say 


The  test  is  truncated  at  n=  A  and  converted  to  a  FSS  test  if  it  has  not  terminated 
before.  Again  using  the  monotonicity  of  the  likelihood  ratio  function  the  test  can  be 
written  as 


TSPRT: 


if  n^A, 


ifn=A, 


yn 


^  Ain)  ■ 
^  Bin)  m 


otherwise, 
^  A:‘(f) 

<  a;‘(t) 


A;’(i)  say 

a;'(b)  5/0'  ^0 

continue 

say  W, 
say 


(3.17) 


1.  Design  Of  Decision  Processors 

The  thresholds  A  and  B  for  the  SPRT  test  can  be  calculated  using 


A  =  —  and  B  =  j8^.  (3.18) 

a 

From  (3.15),  one  can  see  that  the  resulting  false  alarm  and  miss  probabilities 
ase  Pf,  ^  1/A  =  a  and  P„as,  ^  B  =  Bm- 

To  devise  thresholds  for  the  TSPRT,  we  will  split  the  test  into  two  parts. 
Consider  that  the  test  consists  of  a  SPRT  with  thresholds  A,  ^  and  errors  a^,  and 
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a  FSS  test  with  sample  size  h,  threslmld  r  and  enors  a^,  It  was  shown  in  [Ref.  7] 
that  the  errors  of  the  TSPRT  (3.17)  were  bounded  by 


B  ^  B  +  . 

r^tMpn 


(3.19) 


Therefore,  to  obtain  ^  ^  ^upn  ^  ^  split  ttese  errors  into 

the  sums  of  the  errors  due  to  the  SPRT  and  the  errors  due  to  the  FSS  test.  Specifically, 
we  let 


«  =  ^  (1-P>.  (3.20) 

&  =PA  * 

where  Po  and  pj  are  constants  in  [0,1].  By  setting  =  Pfpi  and  /S^  =  pjfi^,  we  can 
design  the  thresholds  A  and  d  according  to  (3.18).  Similarly  letting  a^=  O-pdoi  and 
=  (1-P])0m,  h  and  f  can  be  obtained  using  (3.5).  Note  that  if  Po=Pi=0  the  TSPRT 
becomes  an  FSS  test,  while  siting  both  parameters  to  unity  converts  the  test  a  pure 
SPRT.  Hence  for  values  in  between,  the  TSPRT  can  be  considered  a  mixture  of  SPRT 
and  FSS  test. 

2.  Mean  and  Variance  Analysis  of  The  TSPRT  Scheme 

It  was  shown  in  [Ref.  7]  that  a  properly  designed  TSPRT  was  superior  to  the 
FSS  test  in  the  average  sample  number  (ASN)  for  a  single  search  state.  It  was  also  found 
that  TSPRT  and  SPRT  performed  very  close.  For  the  remainder  of  this  chapter  we  will 
focus  our  attention  to  the  TSPRT  and  make  an  overall  acquisition  system  analysis  using 
the  flow-graph  technique. 


50 


The  system  model  depicted  in  Section  B  of  Chapter  n  mnains  the  same,  only 
the  decision  processor  is  replaced  with  a  TSPRT  schme.  However,  the  time  delays 
associated  with  each  state  are  now  random  variables  instead  of  constants  as  were  in  the 
FSS  scheme.  Therefore  we  will  use  the  statistical  averages,  namely  E{tJ,  E{t^}.  E{t/} 
and  E{t,'^}  in  place  of  the  constant  time  delays  of  the  FSS  test,  to  obtain  mtacq  ^  <’^tacq 
for  this  system.  Recalling  the  equations  from  Clu^er  n  and  using  the  same  notation  and 
symmetry  conditions,  the  required  parameters  for  the  analysis  are 

£(rj  =^5N,.,.7,  =ASN,,.T^, 

EitJ  =  ASN^yT^  =  ASN,yT^, 

~  ■^^^-2.0‘^c  ~  2\\ 

Bfli)  -  J).^^ 

£((];)  = 

ml)  = 

where  ASN  denotes  the  "average  sample  number",  while  VSN  stands  for  "variance  of 
the  sample  number".  The  delays  associated  with  the  coincidence  detector  are 

E{t, )  =  c, .£(/,),  (3.22) 

Eitl^  =  clE{tf) 

where  c,  is  a  known  constant(see  equation  (2.21)).  The  values  of  ASN,  VSN  and  all  the 
detection  and  false  alarm  probabilities  can  be  found  using  simulation  techniques  described 
in  (Ref.  7]. 
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a.  Optimal  Choices  of  Design  Parameters 


For  the  TSPRT  test  we  followed  a  similar  qiproach  to  that  of  the  FSS 
scheme.  However  in  this  case  we  do  not  have  closed  form  e)q>ressions,  therefore,  we 
have  to  d^nd  solely  on  the  simulation  reailts. 

Specifying  the  desired  false-alarm  and  detection  probabilities  suffices  to 
design  a  TSPRT  processor  as  described  above.  We  will  be  looking  for  the  two  design 
parameters  a  and  /3*  which  minimize  the  variance  of  the  acquisition  time. 

We  again  consider  the  acquisition  scheme  dq>icted  in  example  1  of 
Section  B.  The  decision  processor  is  rq)laced  with  a  TSPRT  scheme  and  a  m-sequence 
with  period  N=  1023  is  used  as  the  PN  sequence.  The  mixture  rates  are  taken  as  po  = 
Pj  =  0.5.  The  simulations  are  executed  for  various  a,  jS  pairs,  where  a  6  {0.001 ,  0.01 , 
0.05,  0.1,  0.15,  0.2}  and  jS  €  {0.5,  0.6,  0.7,  0.8,  0.9).  We  have  also  considered 
several  coincidence  detector  parameters  A  and  B  €  [0,8]. 

For  all  the  cases,  all  of  the  actual  values  came  out  to  be  slightly 
greater  than  jS’s,  while  all  the  actual  Pji,  values  were  slightly  less  than  a’s,  which 
indicated  that  the  system  setup  and  :q)proximations  worked  weU. 

Among  all  the  cases  investigated,  the  variance  is  minimized  by  the  pair 
(a’=0.05,  j8*=0.6)  with  the  coincidence  drtector  parameters  A=5  and  5=2.  All  the 
relevant  parameters  associated  with  these  design  parameters  are  given  in  Table  3.4. 

The  notable  improvement  of  the  TSPRT  over  the  FSS  scheme  can  be 
clearly  seen  from  these  results.  Recalling  the  best  case  of  FSS  scheme,  the  variance  of 
the  acquisition  time  is  decreased  from  5715.8  T®  to  2461.4  T*,  (approximately  57  % 
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TABLE  3.4  The  results  of  the  TSPRT  scheme. 


A  =  5  B  =2 

SNR  =  -10  dB 

7  =  A  =  0.5 
o’=  0.05  /3‘=0.6 

E(tJ  =  73.4800  T, 

E(t^2}  =  63.1390  T, 

E(tjJ  =  56.1900  T. 

=  6121.0  V 

E(tJ)  =  4751.0  T.' 

E(0  =  3720.2  T,2 

Mtacq  =  72.1909  T 
Aacq  =  2461.4  T* 

Pjj  =  0.7660 

E(t^/)  =  190.40  T, 

E(t,2^)  =  41096  V 

Pj2  =  0.0840 

E(t^')  =  267.15  T, 

E(t^'^)  =  85057  T,* 

II 

o 

8 

oo 

oo 

E(0  =  230.51  T. 

E(tj^^)  =  62609  T,2 

decrease),  while  the  mean  acquisition  time  is  also  improved  by  iq}proximately  38.3  % 
for  the  same  design  parameters.  We  remark  that  this  improvement  is  a  result  of 
decreasing  the  average  sample  size  per  state,  while  keeping  the  actual  probability  of 
detection  and  probability  of  false  alarm  within  reasonable  values.  Note  that  the  actual  Pj 
is  greater  than  and  the  actual  Pf,  is  less  than  a,  which  is  usually  the  case,  because  the 
system  is  designed  on  a  nominal  worst  case  basis. 

The  minimum  /itacq  for  the  design  example  is  70.1626  T  and  the 
corresponding  a^TACQ=  2814,3  7*.  These  values  occur  for  the  (a,j8)  pair  (0.05,0.5).  As 
we  can  see  that  we  still  have  a  different  set  of  parameters  for  minimum  average 
acquisition  time  choices  as  we  had  for  the  FSS  scheme.  But  the  values  are  somewhat 
closer  for  the  TSPRT  so  that  our  ailment  for  selection  of  the  design  parameters  using 
the  Thebycheff’s  Inequality  holds  even  more  confidently  for  the  TSPRT  scheme. 
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As  one  can  see,  we  again  have  a  small  a  value  and  a  moderate  j8  value 
to  minimize  either  the  mean  or  the  variance  of  the  acquisition  time.  The  ASN  and  VSN 
increases  with  the  increasing  and  decreasing  a  (see  Figure  3.12).  Increasing 
excessively  afTects  ASN  and  VSN,  especially  under  Hq  hypothesis.  Since  a  very  large 
portion  of  the  total  states  fall  under  Hg,  this  results  a  large  mean  and  variance.  Therefore, 
the  system  is  optimized  on  Hg  by  keying  a  low  to  reduce  the  effects  of  false  alarms,  and 
keeping  /3  reasonably  moderate,  not  to  increase  the  ASN  and  VSN  excessively. 

To  certify  the  improvement  achieved  using  a  TSPRT  scheme  over  the 
FSS  scheme,  the  simulations  are  rqreated  for  two  other  per-chip  design  SNR  values  (-13 


Figure  3.12  The  variation  of  ASN  and  VSN  with  various  a  and  /S. 
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dB  and  -3  dB)  with  all  the  other  parameters  remained  fixed.  Figures  3.13  and  3.14 
compare  the  mean  and  the  variance  of  the  acquisition  time  with  the  results  givra  in  Table 
3.2  for  the  FSS  scheme. 

MEAN(T) 


•  ISdB  -lOdB  -SdB 

PER-GHIP  SNR  (dB) 

Figure  3.13  The  resulting  mtacx}  for  various  design  SNR  values. 
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Figure  3.14  The  resulting  o^tacq  for  various  design  SNR  values. 

55 


b.  The  Channel  Mismatch  Problem 


We  examined  the  channel  mismatch  for  the  TSPRT  scheme  in  a  similar 
fashion  that  we  did  for  the  FSS  scheme.  The  system  given  in  Table  3.4  is  designed  using 
the  three  separate  per-chip  SNR  values  (-13,  -10  and  -3  dB)  from  Figures  3. 13  and  3. 14. 
The  effective  SNR  of  the  incoming  signal  is  varied  from  -20  dB  and  0  dB  for  all  the  three 
cases  and  the  plots  for  ^itacq  <^tacq  ^  given  in  Figures  3.  IS  and  3.16  respectively. 

We  see  a  very  similar  picture  to  that  of  the  FSS  scheme.  The  degradation 
caused  by  weaker  signals  are  far  greater  than  the  stronger  signals  and  this  degradation 
is  more  severe  for  higher  design  SNR  (SNRJ.  As  a  numerical  example,  for  SNRd=  0.1 
the  mean  value  is  increased  by  a  factor  of  291  by  a  decrease  of  -10  dB  in  effective  SNR, 
while  a  10  dB  increase,  elevates  the  same  parameter  by  only  a  factor  of  3.44.  For  all  of 
the  three  cases,  slightly  stronger  signals  in  vicinity  of  the  SNR,,  actually  yield  a  slightly 
better  performance. 

As  a  final  observation,  we  can  see  that  the  aforementioned  degradation 
for  the  weaker  signals  seem  to  be  slightly  more  severe  for  the  TSPRT  scheme  than  the 
FSS  scheme,  i.e.,  the  curves  are  steeper.  This  behavior  is  consistent,  since  the  TSPRT 
scheme  is  along  the  optimality  boundary,  therefore  more  sensitive  to  the  external  effects. 

From  all  this  discussion,  we  can  once  again  conclude  that  if  the  system 
designer  has  some  uncertainty  about  what  the  received  SNR  would  be,  it  is  better  to 
design  the  system  using  lowest  possible  SNR,,  to  counter  balance  the  channel  mismatch 
problem.  More  will  be  said  about  the  degradation  in  effective  SNR  in  the  next  chapter 
when  we  examine  the  effects  of  fading  on  the  performance. 


LEGEND:  (DESIGN  SNR)  o— 13dB,  >— lOda  •— IJdB 


EFFECTIVE  SNR 

Figure  3.15  The  variation  of  /xtacq  with  the  effective  per-chip  SNR  for  the  TSPRT. 
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Figure  3.16  The  variation  of  o^tacq  with  the  effective  per-chip  SNR  for  the  TSPRT. 


57 


IV.  PERFORMANCE  IN  THE  FADING  MULTIPATH  CHANNEL 


A.  INTRODUCTION 

The  previous  chapters  have  described  the  design  and  performance  of  the  acquisition 
systems  operating  over  the  classical  AWGN  channel. 

In  this  chapter  we  consider  the  problem  of  the  receiver  performance  for  more 
complex  channels,  namely,  channels  having  randomly  time-variant  impulse  responses. 
This  characterization  serves  as  a  model  for  signal  transmission  over  many  radio  channels 
such  as  shortwave  ionospheric  radio  communication  in  the  HF  frequency  band, 
tropospheric  scatter  (beyond-the-horizon)  radio  communications  in  the  UHF  and  SHF 
bands,  and  ionospheric  forward  scatter  in  the  VHF  frequency  band  [Ref.  15].  The  time- 
variant  impulse  responses  of  these  channels  are  a  consequence  of  the  constantly  changing 
physical  characteristics  of  the  media.  As  an  example,  the  ions  in  the  ionospheric  layers 
that  reflects  the  signals  are  always  in  motion  which  ^>pears  to  be  random  to  the  use'  f 
the  channel. 

Although  fading  arises  from  variety  of  reasons,  fading  phenomena  can  often  be 
modeled  as  causing  multipath  distortion.  That  is,  many  fading  conditions  can  be  modeled 
as  causing  several  alternate  transmission  routes  to  arise  between  the  transmitter  and  the 
receiver.  Since  this  multipath  is  a  result  of  the  time  varying  characteristic  of  the  channel, 
the  nature  of  the  multipath  varies  with  time.  That  is,  in  addition  to  the  transmission  path 
which  the  system  designer  considered,  additional  unplanned  transmission  path  varying 
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with  time  and  appearing  as  random  to  the  user  of  the  channel  are  excited.  The  signals 
arriving  at  the  receiver  through  different  paths  interfere  either  constructively  or 
destructively  to  result  in  signal  fading.  Since  this  variations  appear  to  be  random,  fading 
is  best  described  using  statistical  terms. 

Fading  is  a  widely  studied  area  of  the  communication  theory  and  more  detailed 
explanation  about  the  topic  can  be  found  in  existing  literature. 

Our  intent  for  this  chapter  is  to  examine  the  effects  of  fading  on  our  acquisition 
receiver  performance.  Effect  of  fading  on  the  very  same  system  was  previously  studied 
in  [Ref.  10]  for  a  single  search  state  cell  basis  to  examine  the  variation  of  the  average 
test  length,  the  detection  and  the  false  alarm  probabilities  for  the  TSPRT  scheme.  We 
will  extend  this  study  to  the  overall  system  basis  using  the  flow-graph  technique. 

B.  MATHEMATICAL  MODEL 

The  fading  multipath  channel  under  consideration  is  characterized  as  flat  and  slow 
fading  channel.  Flat  in  this  context  means  that  fading  affects  all  frequency  components 
present  in  the  transmitted  signal  in  exactly  the  same  manner  as  opposed  to  bandwidth- 
selective  channels.  The  received  signal  in  such  a  case  can  be  regarded  as  multiplied  by 
a  random  variable  which  accounts  for  the  fading.  The  other  characterization,  slow  means 
that  the  channel  variations  (which  give  rise  to  the  signal  power  fluctuations)  are  slower 
than  the  lowest  frequency  component  in  the  signal.  Therefore,  the  received  signal  can 
again  be  regarded  as  multiplied  by  a  random  variable,  rather  than  by  a  fluctuating  time 
function  which  would  be  the  case  for  fast  fading  chaimels. 
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In  the  multipath  i^roximation,  the  power  in  the  stronger  signal  path  for  which  the 
system  is  intended  to  be  designed,  is  called  the  direct  component,  while  aU  the  total 
power  in  the  remaining  weaker  paths  is  termed  as  the  diffuse  component. 

It  has  been  shown  in  [Ref.  10]  that  for  the  assumed  flat,  slowly  varying  fading 
channel,  the  received  signal  is  effectively  multiplied  by  a  random  coefficient  which 
has  a  Ricean  probability  density  function,  given  by 

/♦(^)  =  2^(l+r)e-'-*’<'^%(2^V^(I+r)),  ^>0 
where  r  the  ratio  of  the  power  in  the  direct  component  and  the  power  in  the  diffused 
component  (2o^),  with  the  constraint  that  2ff^  =  1.  The  parameters  s^,  and  hence 
r  are  dq)endent  on  the  nature  of  the  channel.  Note  that  for  r  =  oo  (or  j^=l  and  o^=0), 
it  corresponds  to  no  fading.  Also  when  r=0  (or  and  o^=0.5)  we  have  Rayleigh 
fading,  of  which  the  pdf  is  given  by 

4(^)  =  2^e-'^,  ^  >0.  (4.2) 

C.  TEST  STATISTICS  DENSITY  FUNCTION 

Recall  from  Chapter  n  that  the  receiver’s  (Figure  2.1)  test  statistics  is  given  by 
equation  (2.6).  Assuming  a  slowly  varying,  flat  fading  multipath  channel  the  resulting 
test  statistics  can  be  found  as  follows.  The  received  signal  through  the  fading  channel  is 
now  given  as 
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r(r)  =  y4o^a(/  +  iAT^)cos{uj+ff)  +  nit)  (43) 

where  ^  is  the  fading  random  variable  with  the  pdf  of  (4.1)  or  (4.2).  Note  that  since  we 
have  assumed  slowly  varying  fading  channel,  ^  does  not  change  very  much  within  one 
data  bit  interval,  therefore  can  be  assumed  constant  during  the  test  length.  Referring  to 
the  receiver  structure  the  inphase  and  the  quadrature  components  depicted  in  equation 
(2.2)  now  become 

X.  =  ^^pT^S„cose  +  N., 

"  y  (4.4) 

where  and  N^„  are  still  defined  by  equation  (2.3).  As  before,  the  test  statistic  for 
determining  the  alignment  is 

y.  -  xL  <■  xl.  M-s 

The  conditional  probability  density  function  of  y„  conditioned  on  is  non-central 
Chi-squared,  and  given  by 


yn  >0 


(4.(») 


where  K  =  iAjl4)TX^  and  =  nTJSfo  /4. 

It  follows  from  the  above  equation  that  to  detennine  the  test  statistic’s  density 
function,  the  ^  dependance  has  to  be  integrated  out.  Namely,  the  test  statistics  density 
function  is 
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(4.7) 


A.O’J  -  I  di,. 


In  the  previous  study  [Ref.  10]  this  integral  was  carried  out  using  the  probability 
density  functions  of  equations  (4.1)  and  (4.6).  The  resulting  pdf  of  new  test  statistics 
under  fading  conditions  is  rq)eated  below. 


/y  W  = 


l+r 


2<r^[(Ur)+X„/2(r^ 


exp 


(4.8) 


1 

This  result  shows  that  the  probability  density  function  of  y„  is  still  non-central  Chi- 


squared.  The  density  function  is  still  the  same  form  of  equation  (2.6)  but  the  parameters 
X„  and  are  replaced  by  new  values, 


Xl  =  X. 


1  +  r’ 


2' 


2  1  +  r  +X„/2ff„ 
1  +  r 


(4.9) 


D.  PERFORMANCE  ANALYSIS  TECHNIQUES 

Since  we  have  the  precise  description  of  the  receiver’s  test  statistic  density  function, 
we  might  perform  our  ongoing  analysis  methods  to  examine  the  performance  of  the  FSS 
and  TSPRT  schemes  under  fading  conditions.  Note  that  for  the  FSS  scheme  we  have 
closed  form  expressions  and  we  might  perform  the  analysis  analytically.  But  for  the 
TSPRT  scheme  we  have  to  depend  on  simulation  results. 
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In  the  previous  study  [Ref.  10],  the  performance  under  fading  ccmditions  is 
examined  using  a  numerical  integration  technique  to  integrate  out  the  dq)endance  from 
the  conditional  density  function.  That  technique  required  a  number  of  simulations, 
specifically  25  simulations  for  a  single  case  (because  the  area  under  pdf  is  divided  up  to 
25  intervals).  In  this  thesis  we  wUl  use  a  different,  more  direct  simulation  technique, 
that  is,  we  will  design  our  receiver  in  the  normal  way  without  the  fading  and  simulate 
the  effects  of  fading  using  the  pdf  in  equation  (4.8). 

It  is  instructive  to  examine  the  specifics  of  the  simulation  technique  to  verify  the 
validity  of  the  results.  In  the  original  simulation  program  the  test  was  based  on  the 
likelihood  ratio  of  (2.8),  so  that  for  the  SPRT  part  of  TSPRT,  i.e.,  n  <  h  (see  3.16), 
the  test  is  conducted  by 


6=ln(B) 


Ho 

<  Kq-Kx  ^  Jj, 
> 

continue 


W(y>5(X„.,/ab  ) 


> 

< 

continue 


c=ln(i)  (4.10) 


where  A  and  ^  is  calculated  using  (3.18)  and  (3. 20). The  test  statistics  y„  was  then 
normalized  by  a^,  i.e. ,  a  new  random  variable  was  introduced,  z„  =  yjo^.  With  this 
transformation,  the  test  was  written  as 
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(4.11) 


com.  ^  ^  com. 


where  the  thresholds  were  now  functions  of  n.  because  depend  on  n  (see 

equations  (2. 12)  and  (2.13)).  All  (KJo^  values  in  equation  (4. 1 1)  were  calculated  using 
the  nominal  worst  case  values  of  equation  (2. 13).  The  probability  density  function  of  the 
new  test  statistics  z„  was  given  as 


(4.U) 


where  \  under  Af,  and  =  under  Hq.  The  simulation  program 

generates  z„  with  this  statistics,  using  the  actual  values  given  by  equations  (2. 10) 

and  (2.12),  then  calculates  the  value  of  the  middle  expression  in  (4.11)  and  compares 
the  resulting  value  with  the  two  thresholds  a  (n)  and  h  (n)  to  resolve  between  the 
hypotheses. 

We  have  seen  that  the  fading  effects  the  statistics  of  y„  in  the  manner  dq)icted  in 
(4.9).  Using  the  transformation  z„  =  yjo^  this  changes  the  probability  density  function 
of  z„  as 


A(0  =  ,*  .  « h  !  Wo’)  (4-13) 

2(<r„/a/ 

where  the  new  parameters  can  be  written  in  terms  of  the  older  ones  as 
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(4.14) 


x/  _  ^  \  r 

o'h  ,  \*r*\  l2a\ 

(JLf  =  - 

Once  r  is  fixed,  we  can  use  equations  (4.13)  and  (4.14)  to  generate  the  test 
statistics,  therefore,  the  effects  of  fading  on  the  TSPRT  scheme  can  be  examined  for 
various  fading  conditions. 

The  FSS  scheme  under  fading  conditions  is  easier  to  examine.  We  simply  use  the 
cumulative  distribution  function  of  y„  corresponding  to  (4.8)  to  calculate  the  detection  and 
false  alarm  probabilities.  The  cumulative  distribution  function  of  y„  can  be  written  as 


W  =  1  -  Q 


1  +  r  +  X„/2a: 


(l*r)(yyj 


1  \/2al 


(4.15) 


where  the  Q  function  is  defined  as  equation  (3.4). 


E.  RECEIVER  PERFORMANCE  WITH  FADING 
1.  Numerical  Results 

We  now  use  the  flow  graph  technique  to  examine  the  variations  of  the 
performance  parameters  of  the  acquisition  time,  Hjacq  o^tacq  for  various  fading 
conditions.  The  analysis  for  TSPRT  and  FSS  schemes  will  follow  the  stq)s  outlined  in 
Chapter  m.  Subsection  C.2.  and  B.2.  respectively.  We  first  present  the  numerical  results 
for  the  best  case  (minimum  variance)  design  examples  of  Cluq)ter  m  for  SNR  =  -10 
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dB.  The  sy^em  for  TSPRT  has  the  parameters  (a=  0.05,  /9=0.6,  A=5,  B=2)  and  the 
FSS  system  has  the  parameters  (a=0.09,  /3=0.82,  A=S,  B=3).  In  these  figures  the 
fading  conditions  is  varied  for  r  €  {0.0, 0.5,  1.0,  2.0,  5.0,  10.0,  20.0,  50.0,  100.0}  for 
both  of  the  systems. 

Figures  4. 1  and  4.2  show  the  variation  of  the  /itacq  ^  <^tacq  for  the  fading 
conditions  respectively.  Referring  to  these  figures,  we  can  say  that  severe  fading,  i.e., 
r  <  5  cause  a  considerable  performance  degradation  for  both  of  the  schemes.  This 
degradation  is  unavoidable  if  additional  counter  measures  were  not  taken  to  reduce  the 
effects  of  fading.  But  we  can  still  see  that  TSPRT  performs  considerably  better  than  FSS 
scheme  even  under  severe  fading  conditions.  Note  also  that  as  r  increases  both  curves 
asymptotically  approach  to  their  non-fading  values,  which  indicates  that  the  system  setup 
and  analysis  technique  works  weU. 


0  0.5  1  1.5  2  2.5  i  3.5  4  4.5  5 

r  (FADING  VARABLE) 


Figure  4.1  The  variation  of  /xtacq  under  fading  conditions. 
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L£GCNO:  o>1SPRT  »FSS 


r  (FADING  VARIABLE) 


r  (FADING  VARIABLE) 


Figure  4.2  The  variation  of  o^tacq  under  fading  conditions. 

2.  Fading  as  a  "Channel  Mismatch"  Problem 

We  have  examined  the  channel  mismatch  problem  for  both  of  the  FSS  and 
TSPRT  acquisition  schemes  in  the  previous  chapter.  We  have  shown  that  incoming 
signals  that  have  different  strength  than  predicted  by  the  system  designer  caused 
performance  degradation.  Signals  coming  through  fading  multipath  channel  can  be 
regarded  as  a  specific  case  of  channel  mismatch. 

Considering  equation  (4.9),  one  can  see  that  under  fading  conditions,  the 
reductions  to  the  effective  SNR  would  be  due  to  a  decrease  in  the  mean  of  Gaussian 
quadrature  components  of  the  resulting  Chi-squared  test  statistics  (y„ )  and  an  increase 
to  the  noise.  The  parameter  is  decreased  by  a  factor  of  [r  /  (1  -I-  r)]  and  is 
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increased  by  a  factor  of  (  1  +  [(X„  /  (1  +  r)] ).  Since  all  our  formulations  are 

based  on  the  ratio  (VO*  prefer  to  work  on  the  equivalent  imnnalized  test  statistics 
z„(see  (4.12)).  With  the  help  of  equation  (4.14),  we  can  see  that  fading  changes  the 
parameters  of  z„  as 


2 

Oz  = 


r 

under  /f , 

2 

2 

Hr 

\,0 

r 

under 

2 

On 

1+r 

1  1  +  under  H. 

1+r 

1  -*  1  +  under 

1+r  ° 


(4.16) 


Note  that,  since  fading  will  effect  our  decision  scheme  differently  under  each 

of  the  two  hypotheses. 

Consider  first  the  case  under  Hq,  Since  it  rq)resents  a  non-synchronization 
condition,  ideally,  i.e.,  complete  random  sequences  were  used,  we  should  only  detect 
noise  after  the  correlation  process  in  the  receiver.  However,  since  we  are  using  a  partial 
correlation  receiver  and  an  m-sequence,  we  treat  the  resulting  small  correlation  value  as 
if  it  were  a  "signal"  under  Hq.  Our  design  assumes  that  =  (X„.o  lo^  )  is  upper 
bounded  by  per-chip  SNR  (see  equation  (2.13)),  which  is  in  this  case  0.1,  and  fading 
further  reduces  this  value.  Therefore,  we  expect  actually  a  slight  improvement  on  the 
performance,  but  the  slight  increase  in  the  noise  variance  counter  effect  the  performance. 
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Therefore,  we  can  conclude  that  fuling  has  no  or  little  effect  on  the  performance  of  the 
acquisition  scheme  under  Hg.  This  conclusion  will  later  be  confirmed  by  numerical 
results. 

Under  Hj,  we  see  a  completely  different  picture.  In  this  case,  we  actually 
need  a  signal  to  detect  the  acquisition  condition,  and  the  higher  this  signal  value  is  the 
better  the  performance.  Unfortunately,  fading  affects  this  value  adversely;  it  is  decreased 
by  a  factor  of  (r  /  (1  +r)).  In  addition  to  that,  the  noise  variance  (4.16)  is  increased  in 
much  greater  pace  than  under  Hq  (because  X„ ,  >  X„,o). 

Our  design  approximates  the  parameter  ( =  (Xn,,  la2 ))  as  nSNRil  - 1 7 1 

(see  equation  (2. 13)).  Considering  5yVR(l-  [7 1 A)^  as  a  constant,  (X„ ,  /a^ )  is  proportional 
to  the  sample  size  value  n.  Therefore  the  noise  variance  (4.16)  is  also  increased  with 
increasing  n  for  a  fixed  r.  We  expect  that  the  effects  of  fading  to  be  severe  for  H,  under 
these  conditions. 

The  noncentrality  parameter  of  the  resulting  test  statistics  is  effectively 
decreased  and  underlying  noise  variance  is  increased  due  to  the  fading  under  Hi.  Since 
the  incoming  signal  plus  noise  appears  to  have  less  signal  component  and  more  noise 
component  than  no  fading  condition,  the  output  of  the  correlator  in  the  receiver  tends  to 
be  more  "Hq  like".  Therefore,  we  expect  a  decrease  in  the  effective  probability  of 
detection. 

Considering  equation  (4.16),  one  can  see  that  the  worse  the  fading  is  the  less 
the  signal  component  and  the  more  the  noise,  so  that  under  the  worst  case  (Rayleigh) 
fading  the  we  will  have  no  signal  component  but  a  noise  which  the  variance  at  its 
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maximum.  Since  the  noise  variance  also  increases  with  increasing  sample  size  for  a  fixed 
r,  the  output  of  the  correlator  tends  more  and  more  to  be  in  the  extreme  regions  (either 
Hg  or  Hi)  rather  than  in  between,  as  the  test  length  increases.  The  average  sample  size 
(ASM)  is  therefore  expected  to  be  decreased  for  worse  fading  conditions,  i.e. ,  the  TSPRT 
will  reach  a  decision  earlier  (not  necessarily  a  correct  one)  because  of  tl»  increased 
variance  of  the  noise  due  to  fading  conditions. 

We  now  present  the  numerical  results  for  TSPRT  to  validify  our  discussion. 
Figures  4.3  and  4.4  depicts  the  variation  of  the  aforementioned  parameters  under  Hi  and 
Hq  for  a=0.05  and  /3  =  0.6,  0.8,  0.9  and  0.99.  We  can  see  that  the  parameters  under 
Hg,  namely  P^,  and  ASNj.o  essentially  remain  unchanged  under  fading,  while  ASNq  5  and 
more  importantly  Pj  behaves  in  the  way  we  expected.  The  main  cause  of  the  performance 
degradation  under  fading  is  the  degradation  in  the  effective  probability  of  detection.  The 
decrease  in  the  ASN0.5  does  not  help,  because  the  states  under  Hi  is  just  a  small  fraction 
of  the  total  states  (4  out  of  2046  in  our  case),  therefore,  the  reduction  contributes  very 
little  to  the  overall  average  sample  size  per  state  and  the  resulting  mean  and  the  variance 
of  the  acquisition  time. 

Since  the  charmel  mismatch  mainly  occurs  for  Hi  and  decreases  the  actual 
probability  of  detection,  one  corrective  ai^roach  to  reduce  the  effects  of  fading  could  be 
taken  by  increasing  design  j8,  while  keeping  a  the  same.  Since  the  actual  P^  will  also 
increase,  the  system,  therefore  could  satisfy  a  reasonable  detection  rate  and  perform 
better  at  severe  fading  conditions.  But  the  increase  in  j3  will  also  result  in  an  increase  in 
the  average  sample  size  per  state,  therefore,  the  performance  under  no  or  less  severe 
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fading  conditions  have  to  be  sacrificed.  To  illustrate  this,  we  raised  the  value  of  our 
design  example  from  0.6  to  0. 8,0.9  and  0.99.  We  rerun  our  simulation  uncter  the  fading 
conditions  represented  by  r  G  {0  1  5  10  20),  for  these  values  of  while  keeping  all  the 
other  parameters  constant.  We  present  the  results  for  /ijACQ^d  o^tacq  in  Figures  (4.5) 
and  (4.6)  respectively. 


0  10  20  30  40  50  60  70  80  90  100 

r  (FADING  VARIABLE) 

Figure  4.3  The  variation  of  actual  Pj  and  actual  Pj^  for  TSPRT  scheme  under  fading. 
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r  (FADING  VARIABLE) 


Figure  4.4  The  variation  of  ASNq.s  and  ASNj  o  for  TSPRT  scheme  under  fading. 


LEGEND  :  BETA(  o=n.6.  +=0.B.  x=0.9,  *=0.99  ) 


Figure  4.5  The  variation  of  mtacq  under  fading  with  various  0. 


x10<  LEGEND  ;  BETA(  o=0.6.  +=0.8.  »=0.9.  .=0.99  ) 


0  2  4  6  8  10  12  14  16  18  20 

r  (FADING  VARIABLE) 

Figure  4.6  The  variation  of  o^tacq  under  fading  with  various 
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V.  CONCLUSIONS 


In  this  thesis,  we  have  studied  noncoherent  PN  code  acquisition  systems  in  terms 
of  the  mean  and  variance  of  the  acquisition  time  in  the  classical  AWGN  channel  and  in 
the  presence  of  slowly  varying  Ricean  fading.  We  considered  a  single  dwell  fixed 
sample  size  (FSS)  test  decision  scheme  and  a  truncated  sequential  probability  ratio  test 
(TSPRT)  decision  scheme.  The  mean  (/xtacq)  the  variance  (o^tacq  )  of  the 
acquisition  time  was  calculated  using  the  flow-graph  technique.  We  also  assumed  a 
majority  logic  verification  scheme  (coincidence  detector)  with  early  termination  feature 
to  reduce  the  effect  of  costly  false  alarms. 

The  design  of  either  FSS  or  TSPRT  scheme  was  based  on  specifying  the  desired 
false  alarm  probability  and  the  desired  detection  probability,  denoted  by  a  and  j8, 
respectively.  We  have  seen  that  there  is  an  optimum  set  of  (a,/3)  combinations  which 
minimizes  either  /Ixacq  or  o^tacq  in  the  interval  [0,1]  X  [0,1].  The  optimum  value  of  a, 
denoted  by  a*  turned  out  to  be  very  small,  because  the  false  alarms  were  costly  and 
majority  of  the  uncertainty  phases  fell  into  the  category  of  non-acquisition  conditions.  On 
the  other  hand,  the  optimal  denoted  by  /3*  took  moderate  values  (i.e.,  between  0.6  ~ 
0.8).  Our  results  suggested  that  it  was  not  correct  to  assume  that  the  mean  or  the 
variance  of  the  acquisition  time  would  be  minimized  with  (1  -  jS)  «  1.0,  so  that  the 
correct  phase  cell  was  detected  in  the  first  sweep  of  the  uncertainty  region.  Selecting  a 
moderate  0  would  result  a  much  lower  sample  size  per  test  than  a  high  li,  and  this  would 
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result  in  reduced  /xtacq  or  o^tacq  oven  though  several  sweq}s  of  the  uncertainty  region 
might  be  required.  The  TSPRT  scheme  has  been  shown  to  perform  much  better  than  the 
FSS  scheme  in  terms  of  these  performance  parameters.  Approximate  improvement  of  57 
%  in  the  o^tacq  and  38.3  %  in  the  Mtacq  were  achieved  by  the  TSPRT  over  the  FSS 
scheme. 

The  effects  of  various  penalty  times  to  our  design  considerations  were  examined 
for  the  FSS  scheme.  We  have  seen  that  small  variations  in  the  penalty  time  did  not 
interfere  much  in  the  selection  of  design  parameters,  but  an  excessive  increase  in  the 
penalty  time  required  a  much  lower  a,  whereas  fi  remained  essentially  the  same.  We  also 
saw  that  in  case  of  uncertain  penalty  time,  it  was  better  to  use  the  worst  case  value  in  the 
design. 

We  have  examined  the  coincidence  detector  parameters  and  the  early  termination 
feature  of  the  coincidence  detector.  We  saw  that  it  paid  to  employ  a  coincidence  detection 
scheme  in  the  system,  for  it  greatly  reduced  the  mean  and  the  variance  of  the  acquisition 
time.  For  the  optimal  selection  of  the  parameters  A  and  B,  we  found  no  simple  solution, 
because  these  parameters  depended  on  various  system  values.  Our  numerical  results  on 
the  best  cases  suggested  that  A  did  not  have  to  be  so  large,  i.e.,  ^  5,  and  once  A  was 
fixed,  optimum  B  turned  out  to  be  around  \AI2\  ,  depending  on  the  other  system 
parameters.  The  early  termination  feature  saved  substantial  time,  especially  if  a  large 
B  value  (close  to  A)  was  selected.  The  improvement  was  greatly  reduced  for  smaller  B 
values. 
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We  have  examined  the  operating  point  characteristics  of  the  systems  around  the 
design  per-chip  SNR,  what  we  called  the  channel  mismatch  problem.  Our  results  showed 
that  the  performance  degradation  caused  by  weaker  signals  was  far  greater  than  the 
degradation  caused  by  stronger  signals  and  this  degradation  was  more  severe  for  higher 
design  SNR  values.  Following  this  observation,  we  have  concluded  that  if  some 
uncertainty  existed  over  the  effective  channel  SNR,  it  was  safer  to  design  the  system 
using  the  lowest  possible  SNR  value  to  counterbalance  the  chaimel  mismatch  problem. 

Finally,  we  have  considered  a  slowly  varying  Ricean  fading  chaimel  for  both  of  the 
schemes.  The  TSPRT  again  performed  better  than  the  FSS  scheme  under  fading 
conditions.  We  have  seen  that  fading  effectively  reduced  the  actual  detection  probability 
under  the  hypothesis  Hj  while  the  parameters  under  the  hypothesis  Hq  remained 
essentially  unaffected.  We  have  also  seen  that  we  could  improve  the  performance  under 
severe  fading  conditions  by  using  a  larger  design  j8. 
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